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Abstract 

In his study of periodic orbits of the 3 body problem, Hih obtained a formula 
relating the characteristic polynomial of the monodromy matrix of a periodic 
orbit and an infinite determinant of the Hessian of the action functional. A 
mathematically correct definition of the Hill determinant and a proof of Hill's 
formula were obtained later by Poincare. We give two multidimensional gen- 
eralizations of Hill's formula: to discrete Lagrangian systems (symplectic twist 
maps) and continuous Lagrangian systems. We discuss additional aspects which 
appear in the presence of symmetries or reversibility. We also study the change 
of the Morse index of a periodic trajectory after the reduction of order in a 
system with symmetries. Applications are given to the problem of stability of 
periodic orbits. 
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1 Introduction 

In 1886, in his study of lunar orbits, Hill [Ij discovered a formula which expresses the 
characteristic polynomial of the monodromy matrix for a second order time periodic 
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differential equation in terms of the determinant of a certain infinite matrix. Here 
is a slightly modified version of this result. Consider Hill's equation 



a{t)x, (1.1) 



where 



oo 
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i{t) = ^ ate 



k=—oo 

is a real 27r-periodic function. Let p and be eigenfunctions of the monodromy 
matrix. Hill showed that 

=detH, (1.2) 

g27r _^ g-27r _ 2 ' ^ ' 

where H is the infinite matrix^ 



j,kez 



and 5jk is the Kronecker symbol. 

Hill computed det H approximately replacing by a 3 x 3 matrix, which gave 
quite a good approximation. He used equation (jl.2p to find the multipliers approx- 
imately. Astronomical tables obtained by this method are well-known. 

Hill's argument was not rigorous because he did not prove convergence for the 
infinite determinant det H. Several years later Poincare [2\ explained an exact mean- 
ing of the Hill determinant and presented a rigorous proof of Hill's formula. Hill's 
result entered textbooks on differential equations, but was almost forgotten by dy- 
namical systems community until the end of the XXth century when an analogue 
of equation (II. 2p appeared for discrete Lagrangian systems in [3] and independently 
in [1]. Here H turned out to be the finite Hessian matrix associated with the action 
functional at the critical point generated by the periodic solution. In [5] (see also [6] 
a general form of Hill's formula was obtained for a periodic solution of an arbitrary 
Lagrangian system on a manifold. In this case H is a properly regularized Hes- 
sian operator of the action functional at the critical point determined by a periodic 
solution. 

Both discrete and continuous versions of Hill's formula give non-trivial informa- 
tion on the dynamical stability of the periodic orbit in terms of its Morse index. 
Recently this connection was investigated by means of symplectic geometry (see, 
for example, [7| and [8]). However, the approach based on the Hill determinant is 
sometimes simpler and provides additional insight to the problem. 

As mentioned, there are two similar but formally different cases: 



^Hill's matrix was slightly different. 
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— Continuous Lagrangian system with configuration manifold M and r-periodic 
Lagrangian ^{x, x, t) on TM x M which is strictly convex in the velocity. 
Solutions of the Lagrangian system will be called trajectories. Then r-periodic 
trajectories 7 are critical points of the action functional 

^(7)= r ^{^{t),j{t),t)dt 

Jo 

on the set of r-periodic curves 7: M — )■ M. 

- Discrete Lagrangian system with Lagrangian L(x,y) on M x M satisfying 
certain non-degeneracy condition. Then periodic trajectories are n-periodic 
sequences x = (xj)igz which are critical points of the action functional on 
M": 

n 

= ^L{xi,Xi^l), Xi+n = Xi. 

i=l 

Usually one case can be reduced to the other, but this reduction may be cum- 
bersome. Hence it makes sense to consider both cases separately. 

Both versions of Hill's formula look similar. Let P be the monodromy matrix 
of the periodic trajectory, h the second variation of the action functional at the 
periodic trajectory, and H the corresponding Hessian operator. Then 

det(P-/) = fT(-l)™/3det-H', (L4) 

where m = dimM and a = ±1 takes care of orientation. The coefficient /3 is a 
positive scaling factor. 

The operator H is self-adjoint in a proper Hilbert space. For continuous systems, 
H is an unbounded operator, so it needs to be regularized. For example, for Hill's 
equation (jl.ip . H is a, Sturm-Liouville operator. 

Another version of Hill's formula, a generalization of (jl.4p . has the form 

p-"'det{P-pI) = a{-l)"'PdetHp, peC, (L5) 

where Hp is the p-Hessian which coincides with the ordinary Hessian for p = 1. It 
is self-adjoint if \p\ = 1. Since P is symplectic, both sides of (jl.Sp are polynomials 
of degree m in p + p^^. 

Hill's formula (jl.4p has many dynamical applications. The first one is the well 
known statement that the Poincare degeneracy of a periodic trajectory (that is, the 
condition that 1 is an eigenvalue of P) is equivalent to the variational degeneracy 
(the condition det H = 0). 

Another application concerns dynamical instability of a periodic trajectory. It is 
based on the observation that the inequality det(P — /) < implies the existence of 
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a real multiplier (that is, an eigenvalue of P) p > 1. Indeed, F{p) = det(P — pi) = 
det{pl — P) — oo as p ^ +00, and so P(l) < implies the existence of a root 
p > 1. Thus 7 has a positive Lyapunov exponent and is exponentially unstable. 

If det H ^ 0, we have signdetiif = (—1)™"^^, where indH is the Morse index of 
the periodic trajectory. Hence if the periodic trajectory is nondegenerate, then by 
(|1.4|) the inequality cr^ — i^^+md^ < Q implies exponential instability in a 'physical' 
(with /3 > 0) system. 

In some cases it is possible to prove that for any \p\ = 1 the Hessian Hp is 
positive definite and therefore the equation det Hp = has no solutions on the unit 
circle. Then we obtain exponential instability, in fact, total hyperbolicity for the 
corresponding periodic trajectory (Propositions 14.5 1 and I7.2p . 

Below we present other dynamical consequences of Hill's formula. 

Note that the connections between dynamical and geometrical properties of pe- 
riodic orbits are not restricted to Hill's formula. We mention here interesting rela- 
tions between stability properties and the structure (index, signature, and so on) of 
a quadratic first integral of the linearized system ([9] and [IQ]). Many interesting 
results follow from the index formula in symplectic geometry ([7| and [8]). Some of 
our results may be regarded as Lagrangian versions of the results of [11] and |12j . 

Hill's formula is potentially most useful for the study of periodic orbits obtained 
by variational methods. Many such orbits were obtained recently in celestial me- 
chanics by minimization of the action functional on appropriate classes of curves, 
see, for example, |13)- |15] . The most famous example is the figure eight orbit, see 
|13j . However, due to rotational and other symmetries, none of these periodic orbits 
are nondegenerate minimum points of the action. 

In applications periodic trajectories are usually degenerate. For example, any 
periodic orbit of an autonomous continuous Lagrangian system is degenerate. In 
this case the variational equation has a r-periodic solution, 7, and a linear first 
integral, the linearization of the energy integral. Another reason for such a degen- 
eracy (now in both discrete and continuous cases) is the presence of a symmetry 
group, preserving the Lagrangian. This degeneracy also gives r-periodic solutions 
and linear integrals for the variational equation. For degenerate periodic trajectory 
equation (|1.4p is useless because both sides vanish. A nondegenerate version of Hill's 
formula can be obtained with the help of the reduction procedure. We consider the 
case when the Lie algebra V of symmetry vector fields for the variational equation is 
commutative, the dimension of the generalized unit eigenspace of P is 2k, where 
k = dimV (no further degeneracy) and a condition, called the non-degeneracy of 
the trajectory modV, holds. The latter condition has a Lagrangian nature rather 
than Hamiltonian. 

The reduced Hill's formula looks similarly, but the corresponding monodromy 
and Hessian operators P and H-^ act on smaller (reduced) spaces, and are nonde- 
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generate if all the symmetries are taken into account: 

det(P - /) = (7^(-l)™-V^ det H^. 

Here € {1, —1} is the 'reduced orientation' and > 0. 

Now an interesting question appears on the relation between a and a-^ as well 
as between indH and indi?-*-, the Morse indices of the Hessians in the original 
and reduced systems. Indeed, a and indH arc often known for solutions obtained 
by variational methods, while and indi^^ appear in stability problems. The 
following construction explains our answer to this question. 

Let h and be bilinear second variation forms corresponding to the operators 
H and respectively. The forms h and h-^ are defined on the vector spaces X and 
of variations along the periodic orbit, for the original and reduced Lagrangian 
system, respectively. The procedure of the order reduction gives a canonical projec- 
tion U: X ^ X-L. 

For any C ^ ^ r] G X we have h{(^, rj) = 0. Therefore h defines a bilinear 
form hon X = X/V and ind h = ind h. The spaces X and X-^ admit the expansions 

X = e y° e x-^ = n-^®Y-^ 

with the following properties: 

1) dimQ = dimfi-*- = k; 

2) the spaces fi, and Z are /i-orthogonal, while ^l-^ and Y-^ are /i-'--orthogonal; 

3) the restriction is nondegenerate in the discrete case and positive definite 
in the continuous case; 

4) U{Z ) = 0, while the restrictions : ^ Y-^ andU\^: n ^ ^-^ are linear 
isomorphisms; 

5) the forms h\^o ai^d /i"'"|y± coincide in the sense that /i^|yx = ^"""lyx, where 
h = h^ oU; 

6) ^^|n± ~ h^\Q± = X) where x is positive definite in the continuous case. 

In a convenient basis we give an explicit expression for the matrices ^^|qx 
and /i"*" I • Using these expressions we show that 

where the quadratic form b on the generalized eigenspace N = Ker(P— 7)^ is defined 
by b{v) = oj{{P — I)v,v), where w is the symplectic structure. 
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In some cases (—1)™'^^ has a clear dynamical meaning. For example, sup- 
pose that the degeneracy appears solely because the continuous Lagrangian is au- 
tonomous. Then dimV = 1. The periodic trajectory 7 belongs to a smooth family 
of periodic trajectories. Let E and r be the energy and the period along this family. 
Then (see Lemma [6^ 

(-ir^ = -sign^. 

Suppose that the periodic trajectory 7 of an autonomous Lagrangian system has 
no other degeneracy. Then k = 1 and by the reduced Hill's formula it has a real 
multiplier p > 1 provided that 

a(-ir+*'^^^^ <0. (L6) 
ar 

The sign of the quantity dE/dr is easily computed, for example, in the problem 
of the motion of a point in in a homogeneous potential force field. 

It turns out that closed geodesies do not satisfy the condition of non-degeneracy 
mod V . However we show that inequality ()1.6p still implies the existence of a mul- 
tiplier p> I provided no extra degeneracy takes place (Corollary 15. 2p . (Note that 
in this case dE/dr < 0). 

As mentioned above, the subject of this paper is closely related to the theory of 
Maslov-Morse index for periodic orbits of Hamiltonian systems, see, for instance, 
[7], [1], [11], [E]). Some of our results can be obtained by these purely symplectic 
methods. Others are Lagrangian, and so do not have direct symplectic formulation. 
The situation is similar to the relation between Hamiltonian and Lagrangian sys- 
tems: Hamiltonian theory is simpler, more general, and more powerful. Nevertheless 
for many problems the Lagrangian approach is essential. 

This paper splits in two parts: discrete and continuous. Although the majority 
of constructions and statements in the discrete and continuous parts are analogous, 
there are many technical differences which forced us to deal with these two cases 
separately. 

The plan of the paper is as follows. In §[2] we first recall the definition and basic 
properties of discrete Lagrangian systems (DLS). This material is well known to 
specialists, but these objects are not as standard as their continuous analogues. 

Then we present several versions of Hill's formula for a periodic trajectory of 
a DLS. As an application, we give some sufficient conditions for the instability of 
periodic trajectories. Several statements concern stability problem for billiard sys- 
tems in arbitrary dimension. For example, any n-periodic trajectory x of a billiard 
system inside a hypersurface in R'""'"^ such that (_i)™-+"+indx ^ q exponentially 
unstable (Corollarv 12. 6p . As far as we know, there are very few publications about 
stability of periodic trajectories in multidimensional DLS. Here we mention |16] and 
[17] . where trajectories of period 2 are studied. 
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In §[3] we consider DLS with symmetry. We present a discrete version of Routh's 
procedure of order reduction and a reduced version of Hill's formula where the 
degeneracy which appears due to symmetry is removed. We also give a formula for 
the difference between the Morse index of a periodic trajectory of the original system 
and the Morse index of the corresponding periodic orbit of the reduced system. 

In §11] we study reversible DLS, that is, discrete Lagrangian systems with the 
Lagrangian L invariant under time reversal combined with an involution S: M ^ 
M, S"^ = id. Thus L{S{x), S{y)) = L{y,x) for any x,y £ M. For any trajectory 
X = {xi) of the DLS, the sequence x = {Sx^i) is also a trajectory. If x = x modulo 
a translation, the trajectory x is called S-reversible. Then the corresponding space 
of variation splits into a direct sum of spaces of odd and even variations with respect 
to S. Hill's determinant also admits splitting into a product of two determinants. 
Reversible periodic trajectories are also critical points of another action functional 

which is obtained from the original one, by restriction to the space of even 
variations. Morse index of a trajectory with respect to is in general different 
from that computed with respect to . 

Any S-reversible trajectory x of a DLS has 0,1, or 2 fixed points of the involution 
S. According to this we say that typex = 0, 1, or 2. One of application, presented 
in §11] is as follows (Corollarv 14. 4p . Suppose that the billiard surface M C R*""*"^ is 
symmetric relative to a hyperplane and S denotes this symmetry. Let x be an S- 
reversible periodic billiard trajectory of type r € {0, 1, 2} which is a nondegenerate 
minimum of the 'half-length' If m + r is odd, then x is exponentially unstable. 

In §[5]the continuous part of the paper starts. The main technical difference of the 
continuous case is the infinite dimension of the space of variations. Because of this 
the definition of Hill's determinant needs more care. We give a construction defining 
the Hill determinant and present several versions of Hill's formula analogous to the 
ones in the discrete case. Then we give applications to instability of periodic orbits of 
Lagrangian systems including the case of closed geodesies. A typical statement from 
this part (in fact, going back to Poincare) is as follows. Let 7 be a nondegenerate 
closed geodesic on an m-dimensional manifold and a{—\)^~^^^'^'^ > 0. Then 7 is 
exponentially unstable (Corollarv 15. 3p . 

In §[6] we discuss the role of symmetries and give a version of Hill's formula which 
eliminates the corresponding degeneracy. Then we study the relation between the 
Morse index of the periodic trajectory of the original system and the correspond- 
ing periodic solution of the reduced system. We present some applications of this 
formula to the problem of stability for Lagrangian systems with symmetry. 

Finally, in §[7]we consider a reversible CLS. The Lagrangian ^ of an 5-reversible 
NLS is compatible with the involution S in the following sense: 

^{S{x),dS{x)d:,t) = ^{x,-x,-t). 

As in the discrete case, the functional corresponding to even variations is defined. 
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The main questions are the relation between the indices of an 5-reversible periodic 
trajectory with respect to si and s^j^ and the relation between the index with respect 
to and stability properties. We show that in many cases the computation of 
ind 7 mod 2 may be performed on variations from a 2m-dimensional space. 
The authors are grateful to V. V. Kozlov for very useful discussions. 



2 Discrete case 

2.1 Discrete Lagrangian systems (DLS) 

Let M be an m-dimensional manifold and L a smooth^ function on = M x M. 
Denote 

8,L(...) = ^. = ^ (2.1) 

and let 

B{x,y) = -did2L{x,y). 

In local coordinates, 

Bix^y) = -(i^)- (2-2) 



dyj dxi 

In invariant terms, B(x,y) is a linear operator T^M — > T*M, or a bilinear form 
on TxM X TyM. We say that L is a discrete Lagrangian if it satisfies the following 
condition. 

Twist condition. B{x,y) is nondegenerate for all x,y €z M. 
Any discrete Lagrangian L locally defines a map 

T:M^^M^, T{x,y) = {y,z), 

where z = z{x,y) is determined by the equation 

^{L{x,y)+L{y,z)) =d2L{x,y) + diL{y,z) = 0. (2.3) 

In general, T is a multivalued map (relation) with the graph 

r = { (x, y, y, z)GM^xM^: d2L{x, y) + 5iL(y, z) = O} . 

The dynamical system determined by T is called the discrete Lagrangian system 
(DLS) with configuration space M and Lagrangian L. 



Actually, C is enough. 
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Remark 2.1. In this paper we deal with a small neighbourhood of a periodic orbit. 
Hence it is sufficient to assume that the non-degeneracy condition holds locally. 

It is easy to check (see, for example, [18]) that T is symplectic with respect to 
the symplectic 2-form oj = B{x, y) dx A dy, 

uj{u,v) = {B{x,y)ui,V2) - {B{x,y)vi,U2), u = (ni,n2), v = {vi,V2) (2.4) 
((•,•) is the canonical pairing of a covector on a vector). 

Remark 2.2. Let us pass to Hamiltonian variables by the map S: — ?> T*M, 
{x,y) 1-^ {x,px), Px = —diL{x,y). It is locally invertible and replaces T by a locally 
defined map F = STS~^ : T*M — )• T*M. The map F is symplectic with respect to 
the standard symplectic form dpx A dx on T*M, and L is the generating function of 
F: 

F{x,px) = {y,Py), Px = -diL{x,y), py = d2L{x,y). 

Such a symplectic map F is usually called a twist map. 

The map T remains the same after multiplication of the Lagrangian by a con- 
stant, after addition of a constant to L, and after the so-called gauge transformation 

L{x, y) ^ L{x, y) + f{x) - f{y) 

with an arbitrary smooth function / on M. 

A typical example of DLS is the multidimensional standard map: 

L{x,y) = ]^{B{x-y),x-y)-]^{V{x) + V{y)), x,yGM^ (2.5) 

where B \s a. symmetric constant nondegenerate matrix^ 

Consider a domain in R™+^ bounded by a smooth convex hypersurface M. The 
billiard system is a DLS with the Lagrangian L(x,y) = \x — y\ on M x M. Let 
{B{x,y)v,w) be the bilinear form on TxM x TyM corresponding to the operator 
B{x,y) : TxM T*M. A computation gives 

/„/ N \ {v,w) - {v,e){w,e) x-y 

{B{x, y)v, w) = . . , e = r . (2.6) 

\x — y\ \x — y\ 

We may identify TxM and TyM by an isomorphism n(x, y) : TxM — )• TyM, which 
is the parallel projection in ]R™+^ along the segment [x,y]: Ilv = v (mode). Then 

I |2 / \2 

{Bix, y)v, n(x, y)v) = . ~ ^""f ^ > 0, v e TxM \ {0}. 

\x - y\ 

We orient M as the boundary. Since n(x, y) changes orientation, we obtain 



^One can replace the potential {V{x) + V{y))/2 by V{x) or V{y) because they are all gauge- 
equivalent. 
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Proposition 2.1. det B(x,y) < 0. 

Since the image and tlie range of B(x,y) are different, det B{x,y) is not in- 
variantly defined, but its sign is. Tlie fact that the map B(x,y) is nondegenerate, 
provided the hyperplanes T^M and TyM are not parallel to each other in M™', is 
well-known; for a recent reference see [19j. 

In ^18j the reader can find many examples of (mostly integrable) DLS, including 
multidimensional ones. 

For a continuous Lagrangian system (CLS) with Lagrangian ^{x,x), an ana- 
logue of the operator B{x,y) is the matrix ^xxix,x) of second partial derivatives. 
Indeed, consider a DLS on M"^ with the Lagrangian L{x,y) = ^(x, (y — x)/e). In 
the limit as e — > 0, orbits of DLS converge to orbits of the CLS with the Lagrangian 

A computation shows that 

e'B(x,y) = ^ii{x, (y - x)/e) + 0{e). 

In particular, for an analogue of a positive definite Lagrangian system, det B[x^ y) > 
0. 

For m >2 there is no universally accepted discrete analogue of positive definite 
continuous Lagrangian systems. Indeed, in general B is not symmetric and, more- 
over, its symmetry does not have an invariant meaning since B and B* are defined 
on different spaces. Note that the 1-dimensional Aubry-Mather theory was devel- 
oped for twist maps, while multidimensional theory is well developed for continuous 
positive definite Lagrangian systems. 

The most common definition of a positive definite DLS is as follows. Let M = 
and suppose L satisfies the following conditions (see, for example, pQ]): 

- the function L{x, x + v) is periodic in x € T™ and superlinear in v G M™; 

- for any x G the map y diL{x, y) is a diffeomorphism of M™. 

Then L is a generating function of a globally defined symplectic twist map of 
T"" X W^. Evidently, such L satisfies dei B{x,y) > 0. 

2.2 Discrete Hill determinant 

Let (xj)jg2, Xj+„ = Xj, be an n-periodic trajectory of a DLS, that is, r(xi_i,Xi) = 
(xi,Xi+i) for all i. The periodic orbit is determined by x = (xi, . . . ,Xn) € M", and 
a cyclic permutation of x gives the same orbit. By ()2.3p . 

52L(xi_i,Xj) + 9iL(xj,Xj+i) = 0, i = l,...,n, (2.7) 
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where xq = x„ and xi = Xn+i- Thus, x is a critical point of the action functional 

.<2/{x) = L{xi,X2) + L{X2,X3) H h L{Xn, Xl), X G M". 

The point p = {xi,X2) is a fixed point of the map T": M^. The linear 

approximation to dynamics of T near the periodic trajectory is determined by the 
linear Poincare map P = DT'^{p): W — )> W, W = TpM'^. In local coordinates, P 
becomes the monodromy matrix defined uniquely up to a similarity P i— )• S^^PS. 
Eigenvalues of P are called multipliers of the periodic orbit. They determine dy- 
namical properties of the periodic trajectory in the linear approximation. 

Let 

H - ^'^W 

be the Hessian matrix of £/ at the critical point x. Denote 

Bi = B{xi,Xi+i), Xn+l=Xi. 

Theorem 2.1 (discrete Hill formula). 

det(P - /) = ^ J = a(-ir^detH, (2.8) 

Ui=i det Bi 



f7(x) = sign Jl det = H ^^ 



i=l 



1=1 



-1 



(2.9) 



For 'physical' discrete Lagrangians the geometrical meaning of a is the ori- 
entability: the trajectory x is, in a certain sense, orientable if (t(x) > and non- 
orientable otherwise. For example, this is true if DLS is obtained by discretization 
of a positive definite CLS. By Proposition 12. H for a billiard n-periodic trajectory x 
(t(x) = (— 1)". Therefore in this sense billiard periodic trajectories with odd period 
are non-orientable. Note that a is replaced by (— 1)"*"^ if we replace L by —L. 



2.3 Invariant meaning of Hill's formula 

The left-hand side of ()2.8p obviously does not depend on the choice of local coor- 
dinates in M. However an invariant meaning of the right hand side is a priori not 
clear. Let us explain why it is coordinate independent. Let Ei = T^^M. Then 
Bi = B{xi, XiJ^i) is a linear operator Ei — ^ E*^^, and 

Ai = d22L{xi_i,Xi) + diiL{xi,Xi+i) 

is a symmetric operator Ai'. Ei E^. 
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The Hessian of si at the critical point x G M" is a symmetric bilinear form h 
on X = TxM" = E^y. ■ ■ ■ -A En given by 

n 

/i(u,v) = ^Y^[{AiUi,Vi) - {Bi^iUi^i,v^) - {B*Ui+i,Vi)), (2.10) 
1=1 

where 

U = (lii, . . . V = (Ui, ... Uo = Un, Un+l=Ui. 

The form /i is represented by a symmetric operator H: X — > X*: 
/i(u, v) = (Hu, v), u, vGX, 

where 

(Hu)i = AiUi - Bi^iUi^i - B*Ui+i, i = 1,... ,n. 
Define linear operators A, B : X —?■ X* by 

(Au)i = AiUi, {Bu)i = -B,^iu,-i, (B*u), = -B*Ui+i. (2.11) 

Then 

H = A + B + B*. 

Since the Bi are nondegenerate, B is invertible. If we introduce local coordinates, 
then B becomes an (mn x mn)-matrix, and 



detB = (-l)"']^detSi. (2.12) 



i=l 

Hence Hill's formula takes the invariant form0 

det(P - /) = = det(B-^H). (2.13) 

det B 

The equation Hu = gives the variational system of the periodic trajectory x: 

AiUi - Bi^iUi-i - B*Ui+i = 0, G Ei, i e Z. (2-14) 

This is the linear approximation to the system ()2.7p near the periodic trajectory x. 
More precisely, if Ui is any solution of the variational system, then the linearized 
map Pi = dT{xi-i, Xi) acts as 

Pi{ui-i,Ui) = {ui,Ui+i), Ui+i = {B*y^{AiUi - Bi^iUi-^i). 

The kernel of H: X — )■ X* is the set of n-periodic solutions v = (vi), Vi-\-n = Vi, of 
(l2TiD . 

The variational system is a linear Lagrangian system. 



^Since det H, det B are linear operators of 1-dimensional spaces A'""X — > A'""X*, their quotient 
is a well defined scalar. 
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Definition 2.1. A linear periodic discrete Lagrangian system {E,A) is defined by 
n-periodic sequences E = [Ei)i^z of vector spaces and linear operators Ai: Ei ^ E* , 
Bi'. Ei ^ E*j^^, where Ai is symmetric and Bi is nondegenerate. The Lagrangian is 

Ki{ui,Ui+i) = ^ {AiUi,Ui) - {BiUi,Ui+i). (2.15) 
Trajectories of {E,A) are sequences u = (uj) such that 

Thus trajectories of {E, A) satisfy the variational system (j2.14p and are extremals 
of the quadratic action functional 

1 " 

-/i(u,u) = ^Ai(ui,Ui+i). (2.16) 

i=l 

The system {E, A) is the linearization of (M, L) at x. 
2.4 Generalized Hill determinant 

Let us define a generalization of the Hessian H. Let 5^ = {p€C:|p| = l}. 
For any p £ S^, let Xp be the space of all quasiperiodic complex vector sequences 
u = {uj)j£z such that iij+„ = puj. Here Uj lies in the complexification of Ej which 
we will denote Ej for simplicity. The Hessian of the action defines a Hermitian form 
on Xp-. 

n 

v) = - {Bj_iUj_i,Vj) - {B*Uj+i,Vj)). 

i=i 

Since a quasiperiodic sequence {ui)i^i is determined by u = (Mi,...,n„) G 
£^1 X • • • X En = X, we identify Xp with X (more precisely, with the complexification 
of X). Then we obtain the Hermitian form 

/ip (u, v) = (HpU,v), u,vGX, 

where Hp : X — > X* is given by 

(Hpu)j = AjUj - Bj_iUj-i - B*Uj+i, 

j = l,...,n, Uo = p~^Un, Un+l=pUi. 

Similarly we define an operator B^: X X*: 

(Bpu)j- = -Bj^iUj^i, j = 1, . . . , n, uo = p~^Un. 

Then 

Hp = A + Bp + b;. 

The operators Hp, Bp make sense for any non-zero p € C. 



14 



Theorem 2.2 (generalized Hill formula). For any p G C 



det(P - pi) 



det Hp 
detBp 



(2.17) 



Since det Bp = p 



detB, we obtain 



p-"" det{P-pI) 



det Hp 
detB ■ 



(2.18) 



This is an invariant version of the result of [3] . For p = 1, (|2.18p gives (|2.8p . 

Both sides in ()2.18p are polynomials of degree m in the Hill discriminant p + p"^ 
with senior coefficient 1. Indeed, the characteristic polynomial F[p) = det(P — pi) 
of the symplectic operator P satisfies F{p) = p'^^^F{p~^). Hence p~"^F{p) is a 
symmetric polynomial in p and p~^. Thus it is a function of p + p~^. 

In coordinates, Hp is an (mn x mn)-matrix which coincides with H with two 
exceptions: in the upper right (m x m)-block, —Bn is replaced by —p~^Bn and in 
the lower left (m x m)-block, —B* is replaced by —pB*. 

Proof Theorem \2.l^ Let us show that G{p) = det(B~^Hp) is a polynomial of de- 
gree 2m with senior coefficient equal to 1. 

Let = and make a change of variables u i— )■ w, Uj = u^wj. Then the 
operator Hp is replaced by Hpi X — > X*, where 

(lipW)j = AjWj - ly^^Bj^lWj-i - l^BjWj+i, Wn = Wq, Wn+l = Wi. 



Similarly, Bp is replaced by Bp = B. By the invariance of a determinant. 



is a polynomial of order 2mn in u. Since Bp and Hp are linear in p and p ^, G{p) 
is a polynomial in p, p~^. Thus, G{p) is a polynomial of order 2m in p. The senior 
coefficient is det(B-iB*) = 1. 

We have HpU = if and only if u satisfies the variational system and Un = puQ, 
Un+i = pui, or Pw = pw, where w = {uo,ui). Hence F{p) = det(P — pi) = is 
equivalent to det Hp = 0. Thus the polynomials F{p) and G{p) have the same roots 
and so they coincide. □ 



Hence 



H. = A + z^-^B + i/B*. 



G{p) = det(Bp ^Hp) = det(i/B"^A + / + i^^B'^B*) 
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2.5 Some applications 

Identity (|2.8|) implies that dynamical non-degeneracy of a periodic trajectory is 
equivalent to the geometric non-degeneracy: det H = <^ det(P— /) = 0. Actually, 
the proof of Theorem 12.21 was based on this fact. 
Equation (|2.8|) gives 

cj(x)(-irdetHdet(P-/) > 0. 

Corollary 2.1. Suppose that a{x){-l)'^+''''^'^ < 0. Then X is dynamically unstable: 
there is a real multiplier p > 1. 

For example, the hypothesis holds when (t(x)(— I)*" < and x is a nondegenerate 
local minimum of the action £/. 

Corollary 2.2. Suppose that (t(x)(—1 < and x is a nondegenerate local 
maximum of the action £/. Then x has a real multiplier p > 1. 

Indeed, it is sufficient to use the following 

Proposition 2.2. If det{P — /) < 0, there is a real positive multiplier p > 1. 

Proof. Consider the characteristic polynomial F{p) = det(P — pi). Its roots are 
the multipliers of the periodic solution x. We have F(+oo) = -|-oo and F{1) = 
det(P — /) < 0. Then there exists a real root p > 1. □ 

Corollary 2.3. // cr(x)(-l)'°^"-i < 0, there is a real positive multiplier p < —1. 

Indeed, for p = —1 Hill's formula (|2.18|) gives 

cr(x) det H_i det(I + P) > 0. 

Let x^ be the iterate of a periodic trajectory x, that is, the corresponding 2n- 
periodic trajectory. 

Corollary 2.4. Suppose fj(x)(— < q. Then x is exponentially 
unstable. 

Proof. Since 2n-periodic vector fields along x^ are split into n-periodic and n- 
antiperiodic ones, indH(x^) = indH(x) + indH_i(x). It remains to use Corollary 
[Ql □ 

In the case m = 1 there is a possibility to identify hyperbolicity or ellipticity of 
a periodic trajectory in terms of the index. 

Corollary 2.5. Suppose that m = 1. Then a nondegenerate periodic trajectory x is 
hyperbolic if and only if ind x^ if even and elliptic if and only if ind x^ is odd. 
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Proof. The hyperbolicity of x is equivalent to the hyperboHcity of x^. For m = 1, 
x^ is hyperboHc if and only if it has a multiplier p > 1. This is equivalent to the 
inequality cr(x2)(-l)^+™'^^(''^) < 0. It remains to note that cr(x2) = cr^(x) > 0. □ 

Consider the convex billiard bounded by a hypersurface M in M"^^. Then the 
corresponding action is length and, by Proposition 12. H deti3(a;,y) > 0. Therefore, 
a = (—1)" and we obtain 

Corollary 2.6. Suppose (_i)'"+»^+indH(x) ^ g_ xhen x is exponentially unstable 
by Corollaru \2.1\ 

In particular, x is exponentially unstable in each of the following two cases 

• if 771 is odd and x is a nondegenerate local maximum of the billiard length 
functional; 

• if m + n is odd and x is a nondegenerate local minimum of the billiard length 
functional. 

For m = 1 by the Birkhoff theorem [2_lJ (see also [6j), any convex billiard system 
has (at least) two periodic trajectories of period n with rotation number k < n, 
where one of them is a maximum of length, and hence generically hyperbolic. The 
other has index 1, and so det(P — /) > 0. This implies that the trajectory has no 
real multipliers > 1. Indeed, if such a multiplier exists, then the other one is also 
real and greater than 1. This contradicts detP = 1. 

The problem of stability for billiard trajectories of period 2 is systematically 
studied in the recent paper \T7\. In this case, the characteristic polynomial, as a 
function of p + p~^, can be presented as a determinant of some m xm matrix. This 
matrix is explicitly determined by the matrices of second fundamental forms of the 
surface M at the end points of the trajectory. 

The requirement for n to be even in the hyperbolicity condition for a periodic 
trajectory of minimal length (m = 1) at the first glance looks somewhat strange 
because the billiard trajectory minimizing is naturally associated with a locally 
shortest closed geodesic on a two-dimensional Riemannian manifold. Such geodesies 
due to Poincare f2] are known to be hyperbolic. However one should keep in mind 
that this Poincare's result is valid only for orientable geodesies (see details in §[5]) 
while a periodic billiard trajectory with an odd period should be associated with a 
non-orient able geodesic. 

A simple example of an elliptic action minimizing billiard trajectory with odd 
period can be constructed as follows. Let the billiard curve be an acute-angled 
triangle ABC. Then by a well-known theorem from planimetry the projections A', 
B', and C of the vertices to the opposite sides form a triangle (the orthotriangle) 
which presents a local nondegenerate minimum of the billiard action (Fig. [H a). 
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The corresponding periodic trajectory is parabolic: its multipliers pi^2 are equal 
to -1. 




a) b) c) 



Figure 1: Deformation of a parabolic periodic orbit 

A small deformation of the billiard curve does not destroy the periodic trajectory 
A'B'C and just slightly deforms it. If the boundary curve becomes concave, we 
obtain a Sinai billiard [2^2j. In this case the trajectory is hyperbolic (Fig. [H b)). If 
the boundary curve becomes strictly convex (the curvature gets positive (Fig. HJ c)), 
then the trajectory becomes elliptic still having a locally minimal action provided 
the deformation is small. 

3 Continuous symmetry in a DLS 
3.1 Discrete symmetry 

A diffeomorphism tp: M ^ M is a. discrete symmetry of the Lagrangian L if the 

map Tp = Tp X tp: — > M^, ip!{x,y) = (^'ip{x),ip{y)^ preserves L: 

L{'ip{x),ip{y)) = L{x,y). 

A more general definition is that tp preserves L up to a cocycle: 

L{^P{x),^P{y)) = L{x, y) + f{x) - f{y). 

If there exists a function g such that g o -p — g = /, then -p preserves the gauge- 
equivalent Lagrangian L{x,y) = L{x,y) + g{x) — g{y). 

Proposition 3.1. A symmetry takes a trajectory of a DLS into a trajectory. Thus, 
'p o T = T o tp. 

Proof. Since tp preserves the action functional j^, it takes critical points to critical 
points. □ 
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If a DLS (M, L) admits a discrete symmetry group F, then, in principle, sym- 
metry can be removed by a factorization M = M/T of the configuration space. 
However it is useful to keep in mind the following two aspects. 

1. Since in general L{x,y) / L{gix, g2y) for gi ^ g2, the Lagrangian becomes 
multivalued after the factorization M/T. This phenomenon is effectively used in 
the construction of a symbolic dynamics by the method of anti-integrable limit, see 
[25] and [M] (a more general setup is discussed in [23], a continuous analogue is 
presented in [26]). 

2. A periodic trajectory of the original system can turn into a trajectory of the 
factorized system with a smaller period. Therefore the trajectory can lose orientabil- 
ity. Moreover, the configuration space itself can lose orientability. This happens, for 
example, in the case of a billiard system inside a convex hypersurface M C M™"^^ 
symmetric with respect to the origin. Then G = {id, S}, where S{x) = —x. Then 
M = M/G is homeomorphic to the m-dimensional projective space which is non- 
orientable for m even. 

3.2 Noether symmetry 

Let w be a smooth vector field on the configuration space M and tps'- M ^ M its 
phase flow. We say that w is a symmetry field for the DLS if ij^s is a symmetry for L 
for all s. 

Define the vector fields wi = (w, 0) and W2 = (0, w) on M^. Let w = wi -|-W2 = 
(w, w) be the vector field corresponding to the group action tpg = ipg x tpg- We have 
an equivalent version of the definition: w is a symmetry field for L if and only if 

D^L{x,y) = D^f{x)-D^f{y) (3.1) 

for some function / on M. If we replace L by its proper calibration, equation (j3.ip 
can be replaced by 

D^L = D^^L + D^^L = 0. (3.2) 
Proposition 3.2. Let w be a symmetry field for L . Then 

f = D^,L = -D^^L 
is a first integral of the corresponding DLS, that is, ^ oT = ^ . 
Proof. Suppose that {y,z) = T{x,y). Then 

c/{x,y) - /{y,z) = D^^L{y,z) - D^^L{x,y) 

= D^iHy^ + D^iHy, z) = D^L{y, z) = 0. 
The last expression vanishes by (|3.2p . □ 
We call ^ the Noether integral. 
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Proposition 3.3. Let w be a symmetry field. Then the group action ips '■ — )• 
preserves ^ . Equivalently, w is tangent to the level surfaces 

Nc = f-\c) C M^. 

Proof. The derivative of ^ along w is 

D^/ = (Dwi + D^,)D^,L = 

by (j3.2p because tlie differential operators D^^ and Dwi + commute. □ 

Note that ^ is the Hamiltonian generating the group of symplectic trans- 
formations with respect to the symplectic form (j2.4p on M^. 

3.3 Routh reduction of order 

Suppose system (M, L) admits commuting independent symmetry fields w^, . . . , w*^: 

[w",w^]=0, a,P = l,...,k. 

Then the flows "0° : M ^ M of symmetry fields commute. Let G be the correspond- 
ing commutative group acting on M by 

In general the flows "0^^ may be incomplete, and then G is a local group acting on 
M. Since we are interested in a neighbourhoo^of a periodic orbit, these non-local 
questions are irrelevant for us. Suppose that M = M/G is a smooth manifold and 
tt: M — 7- M a smooth fibration (at least locally this is always true). 

Let be the Noether integral of (M, L) corresponding to w" and let ^ = 
(^^, . . . , be the corresponding vector integral. We fix the value c € M'^ and 

restrict T to the level set Nc = ^~^{c) C M^. By Proposition 13.31 the group G 
acts on Nc- If M is a smooth manifold, then A'^ = Nc/G is also a smooth manifold, 
and T defines a map T: N ^ N. It is symplectic with respect to the quotient 
symplectic structure uj on N . We would like to represent T as a discrete jLagrangian 
system with quotient configuration space M and Lagrangian L on M x M. For this 
reduction we need 

Non-degeneracy assumption. The matrix G = {g"^^), 

g'^^ix, y) = {B{x, y)w"(x), w'^(y)), x,y e M, (3.3) 
is nondegenerate. 
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First let c = 0. Let f{x,y) G be the critical point of the function s i— )• 
L{x,ips{y)), provided it exists and is unique. Note that the Hessian of this function 
equals G{x,Tps{y)), and so is nondegenerate. The reduced Lagrangian is defined by 

Z(x,y) = L(x,V'/{a;,y)(y))- (3.4) 

Since L{x,y) = L{x,y) depends only on 5; = 7r(x) and y = 7r(y), it is a function on 
M X M. 

Suppose now that c = (c^, . . . , c^) ^ 0. Locally there exist smooth functions (jjp, 
/3 = 1, . . . , /c, on M such that D^acj)^ = S'p. In general there are topological obstruc- 
tions to the existence of single valued globally defined (j)a- However, if vr: M ^ M 
has fibre M'^, then (p^ exist globally. Since we work in a neighbourhood of a periodic 
orbit, this is irrelevant for us. 

Replace L by gauge-equivalent Lagrangian^ 

L{x,y) = L{x,y) + c'^{<f>a{y) - (f>a{x)) ■ 
Then is replaced by 



Now c = and so L can be defined by ()3.4p with L replaced by L. 
Here is Routh' Theorem for discrete Lagrangian systems. 

Proposition 3.4. The projection n: M ^ M takes trajectories of the system (M, L) 
with ^ = c to trajectories of the reduced Lagrangian system {M,L). 

Next we give a coordinate version of the Routh reduction. Since the result is 
local, it is sufficient to perform the reduction near a given trajectory x'' = (x^). Since 
the vector fields w" are independent and commute, in a neighbourhood Ui of the 
point x^ there are local coordinates yi € M'"^'^, Zi € M*^ such that w°|t/. = d/dzf. 
Similarly to the continuous case coordinates zf are called cyclic. The variables 
yi G W^~^ are local coordinates on M. Equation p.2p means that 



L{xi,Xi+i) = ^{yi,yi+i,Ui), Ui = Zi+i - Zi. 
By p.Sp . the matrix 



^We use Einstein's sum rule with respect to repeated Greek indices, but not Latin indices. 
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is nondegenerate. 

Without loss of generality we assume that c = 0. Then 

/ =du^^{y^,yi+uui) = {). (3.6) 

Equation (|3.6p can be locally solved with respect to Ui = fi{yi,yi+i). Then the 
Routh function is defined by 

L{yi,yi+i) = ^{yi,yi+i, fiiyi,yi+i)). (3.7) 

Hence 

dyl{yi,yi+i) = dy^^{yi,yi+i,Ui)\^_j,^y^^y^^^y (3.8) 

and similarly for the derivative with respect to yi+i- 

Suppose Xi = {yi, Zi) is a trajectory of the system (M, L) with ^ = 0. Then 

dy,{-^{yi-uyuui^i) + ^iyi,yi+i,ui)) = 0. 

By mi, 



dyi {L{yi-i,yi) + L{yi, yi+i)) = o, 

so (yi) is a trajectory of the reduced system (M, L ). 

To finish the proof of Routh's Theorem, it remains to show that L is a discrete 
Lagrangian, that is, it satisfies the twist condition. This follows from 

Lemma 3.1. Let Bi = B{yi,yi+i) = -di2L{yi,yi+i) and Bi = B{xi,Xi+i). Then 

Proof. Putting (ga/Bi) = G^^ we have 

d 52 
— /„(2/„y,+i) = ^-^Jf(y,,2/.+i,tx,)L^^^^(^^ ,^^^^). (3.10) 

We differentiate (j3.8p with respect to yi+i using (I3.10p : 

-Bi = - — h 7^ ^ ^ ^ = 7^ — 7\ h go 



dyi dyi+i %+i duf % % dyi+i " ' %+i dzf_^^ Qy. dz^ ' 

or in the matrix notation Bf' = Bj^ — B"^-^ gafnB^^ . It remains to make an exercise 
in linear algebra. □ 
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3.4 Symplectic reduction for the Poincare map 

Suppose that the periodic trajectory is degenerate. Then the Hnear Poincare map 
P: W ^ W has a unit eigenvalue: there exists w ^ such that Pw = w. Since P is 
symplectic, oj{w,Pu) = uj{w,u), and so Jw{u) = uj{w,u) is a linear first integral of 
P. Then it is possible to reduce P to a linear symplectic map P : W ^ W oi lower 
dimension. 

This section deals with symplectic linear algebra, and the origin of the symplectic 
map P is irrelevant. In particular, the notations below will be used both for discrete 
and continuous Lagrangian systems. 

Suppose there are several eigenvectors corresponding to unit eigenvalue. Let 
V C {w ^ W : Pw = w}. Then P has a first integral J : W ^ V*: for w ^ V , 
{J{u),w) = Jw{u)- We assume that V is isotropic: uj\y = 0. Then V C J^^{0). We 

put W = j~i(o)/y. 

Proposition 3.5 (Poincare). P generates a reduced symplectic operator P: W ^ W 
such that the diagram 

p 

w—^w 

is commutative. Furthermore, 

det(P - plw) = (1 - pf' det(P - pl^), k = dimV. 

3.5 Routh reduction for linear discrete Lagrangain systems 

Next we translate Proposition 13.51 to the language of the variational system, that is, 
the linear Lagrangian system {E,A). To any eigenvector w of the Poincare map P 
there corresponds a non-zero n-periodic solution w = (wi) of the variational system. 
To the periodic solution w there corresponds a linear periodic first integral 

Ij{uj,Uj+i) = {BjWj,Uj+i) - {BjUj,Wj+i). 

Indeed, if u = (uj) is a solution of (|2.14|) . then 

= {AjWj — B*Wjj^i - i?j_iWj„i, Uj) — {AjUj — BjUjj^i - Bj^iuj^i, wj) 

= Ij{Uj,Uj + i) - Ij_i{Uj_i,Uj). 

Hence 

Ij-i{uj_i,Uj) = Ij{uj,Uj+i). (3.11) 
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In fact, Ij{uj,Uj^i) = Jw{u), where is the integral of the Poincare map and 
u £ W corresponds to the trajectory (%). 

Suppose now there are several eigenvectors and let V C Ker(P— /) be an isotropic 
subspace. Denote by F C X the set of periodic trajectories corresponding to V. Let 

G y, a = 1, . . . , /c, be a basis in V. Then the corresponding periodic trajectories 
w" = (wf) form a basis in 7. Let 

Ij'{uj,Uj+i) = {BjWj,Uj+i) - {BjUj,Wj^i) 

be the corresponding integrals of the variational system. Since V is isotropic, the 
integrals commute: 

If{w^,w^+,) = {Bjwf,w^^,) - {Bjw^,wf^,) = u{w-,w^) = 0. (3.12) 

We sometimes write Ij = {ij, . . . , Ij). 

Below we need several non-degeneracy conditions. 

Condition A. The symmetric matrix 

is nondegenerate for all i. Denote (ga/Si) = i^gl 
Condition B. The matrix 

n n 

G = ^Gr^ = (gap), gap = ^g^pi, (3.13) 

4=1 i=l 

is nondegenerate. 

Many of our results hold without condition [Bl so we impose it later. Conditional 
is used almost everywhere, so we impose it now. In the case of CLS, an analogue of 
condition [A] is also introduced, but finally it turns out to be unessential, see ^ IA.3[ 
An analogue of condition [B] is always satisfied for CLS. 

Denote 

Fi = {wi : w G r} = SY>aii{w}, . . . , ) C Ei. 

The reduced Poincare map P corresponds to the reduced linear Lagrangian system 
(£',A) with Ei = Ei/Fi which is obtained by the Routh reduction of the system 
(£", A). Under the non-degeneracy condition [XI dimFj = k and the reduced config- 
uration space Ei = Ei/Fi can be identified with 

E^- = {ue Ei: {Bi_iw°Li,u) =0, a = 1, . . . , /c} 
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via the projection Hi: Ei ^ E^. 

UiU = Ui- gai3i{Bi_iwf_i,u)w^ . (3.14) 
We represent any vector Ui G Ei as 

Ui = Vi + X^iw'^ , Vi = UiUieE:^, X^i = go,j3i-i{Bi_iw1_^,u) . (3.15) 
The Routh reduction for DLS is described by the foUowing 

Theorem 3.1. Let u = {ui), Ui E Ei, be a trajectory of the system {E,A) such 
that Ii[ui,Ui^i) = 0. Then v = (vi), Vi = HjUj G E^ , is a trajectory of the linear 
Lagrangian system {E^,A-^) with the Lagrangian 

A:tivi,Vi+i) = ^ {AiVi,Vi) - {BiVi,Vi+i) - ^ {dvi^Vi), 

where 

{CiVi,Vi) = gai3i{BiVi,wf^^){BiVi,w^^-^^). 

Conversely, ifv is a trajectory of the system (£'"'-, A"*"), then there exists a trajectory 
u of the system {E,A), defined mod F such that Ii{ui,UiJ^i) = and IIu = v. 

For the proof we will need an evident 

Lemma 3.2. Let u, v he such that Ui — Vi (z Fi: 

Ui = Vi + XpiW^. 

Then If{ui, Uj+i) = c"^ for all a = 1, . . . , k and all i if and only if 

AXai = Xai+l - Xai = gafiiic'^ " /f (Uj , Uj+i )) . (3.16) 

Equation (|3.16p follows from 

= If{ui,Ui+i) = Li{vi,Vi+i) +g^^AXfji. 

Proof of Theorem \3.1[ Let u = (uj), Ui G Ei, be a trajectory of {E,A) such that 
Lf{ui,UiJ^i) = 0. Then for any variation (pi G Ei such that = except for 
i = 1, . . . , n, we have 

n 

= /i(u,</)) = '^{AiUi - B*Ui+i - Bi_iUi-i,(j)i). 
1=1 
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Let V = nu. By ([XT5]) and (l3l^ . 

If (fi, Vi+i) = -{BiVi,wf^i), AXai = gal3i{BiVi, wf+i). 

Choose the variation (j) such that (/){ S E^. Using 

A^wf = Bi^iwf_^ + (3.17) 
and {Bi —iWj^_^, (pi) — 0, we obtain 

n n 

h{u, (j)) = ^{AiVi - B*ViJ^i - Bi^iVi^i, 4>i) - ^^i3i{Bi4>i,wf_^_^) 

1=1 1=1 
n 

= ^(v,0) -^gai3i{BiVi,Wi^-^){Bi4)i,w^j^-^) = h^{v,(j)), 

i=l 

where 

n 

/l-L(v,v) = /l(v, v) - ^5'a/3i(^i^'i,W^iVl)(^i^'i,W^f+l)- (3-18) 
1=1 

This is the quadratic action functional for the system {E^,A-^). Since (pi € E;^, 
i = 1, . . . ,n, are arbitrary, v is a trajectory of A-*-). We skip the proof of the 
converse. □ 

For n-periodic v the bihnear form h-^ on X-^ = E^ x • • • x Ej^^ equals 

n 

h-^{-V,y/) =^^{Afvi- Bl*Vi+i- Bl^^Vi^l,Vi), Vn+l=Vi, Vq = Vn, (3.19) 

1=1 

where 

Ai = Ri{A-Ci), Bl = R,+iBi 
and i?i : E^ — )> {E^)* is the restriction map. The reduced variational system is 

Afvi - B^*v,+i - B^_^Vi^i = 0. 

Let H^jB"*-: X-^ — )• X^* be the corresponding linear operators. We also put 

a^P^=(jldetB^^ , /3^>0, a^e{l,-l}. 

Then analogously to (j2.12p . 

detB^ = (-l)'"-^' JJdetS,^. 
Hill's theorem 12.11 applied to the reduced system, gives 
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Corollary 3.1. The following reduced Hill formula holds: 

det(P-/) = = ^^(-ir-'/3^detH^. 

To use this formula for stability problems, we need to know a-^ and the Morse 
index of h'^. However, the relation between the Morse indices of h and /i"*" is not 
evident. The reason is that a periodic sequence v € in general corresponds to 
a non-periodic sequence u ^ X such that Ii{ui,UiJ^i) = 0. We discuss this problem 
in the next two sections. 

Lemma |3. II implies ind-Bj = ind-B,^ + indGj (mod 2). Therefore, 

3.6 ( 

Degeneracy of $h$) Degeneracy of h We denote hy T d X the space of periodic 
solutions corresponding to V C Ker(P — /). It is spanned by wi, . . . ,Wfc € X. Since 
Hw = for w G r, the Hessian bilinear form /i(u, u) = (Hu, u) is degenerate and 
defines a bilinear form h on X = X/V. To compare h with h^ we need to restrict h 
to the level set of /. Let 

y = {ueX : /i(ni,n2) = • • • = J„(u„, txi)}, (3.21) 
Z = {v € X : G Fi} = {v € X : = Kiwf}. (3.22) 

Proposition 3.6. The spaces Y and Z are h-orthogonal, that is, h{u, v) = for 
all u €zY and v € Moreover, Y is the h-orthogonal complement to Z: 

Y = {ueX : /i(u, v) = for all v G Z}. 

The restriction of h to Z is 

n n 

/i(v,v) = Y^gfAX^^AXf,, = Y,{GiAXi,AX,), = X^iwf. 

1=1 i=l 
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Proof. Let u. £ Y and v = (Aajwf ) € Z. Then 

n 

/l(u,v) = ^(^jUj - B*Ui+l - Bi^iUi-l,XaiWf) 

1=1 

n 

1=1 

n 

= ^Xai{{AiWf,Ui) - {BiUi,Wi^-^) - {Bi^iWi_i,Ui)) 

1=1 

n 

= Y,{A^wf - B*wt+i - Bi^iwf_„ X^m) = 
1=1 

by (j3.17p ; we used that 

{BiWi,Ui+i) = {BiUi, + c", {Bi^iUi-i,Wi) = {Bi.iwf_i,Ui) - c". 

Conversely, if v = [Xaiwf) and /i(u, v) = for ah A^, then 
{AiUi - B*Ui+i - Bi^iUi^i,wf) = 0. 

Using (|3.17|) we obtain 
Thus, u G y. 

Next we compute the restriction of h to Z. Let v = (Aajw"). Then by (|3.17p . 

n 

^(v,v) = ^{XaiXi3i{AiWi ,w'^) - XaiXi3i-l{Bi^iwf_^,w'^) 
1=1 

- XaiXpi+i{BiWi ,w'^_^_^)) 

n 

= ^{Xai{Xpi - Xi3i+l){Biwf,W^_^_^) 
1=1 

+ Xpi{Xai - Xi3i-i){Bi_iw'^_^,wf)) 

n 

= ^(^r ^aiiXpi — Xpi^i) + gffiXai{Xi3i — A^j_i)) 
1=1 

n 

i=l 
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□ 

We obtain a quadratic form on Z: 

n 

/i|^(v,v) = (H^v,v) = ^(GiAAi,AAi) = (GA,A), 

where Aj e R*^, : 2" ^ Z*, and the operator G : M'^" ^ M'^" is defined by 

(GA)i = Gi_iAAj_i — GjAAj. 

We have 

Ker G = {A e M'^" : Gi AAi = • • • = GnAAn}. 

n 

Thus, AAj = G^^c. For A € Ker G the equation ^^AAj = imphes Gc = 0, where 

G is the matrix ([333]). _ 

Now we impose the non-degeneracy assumption |Bl the matrix G is nondegener- 
ate. Then 

Ker G = {A G M"*^ : Ai = ■ ■ ■ = A„} and Ker Hz = T. 

Let Z = Z/T. 

Proposition 3.7. The form h^^ ^•^ nondegenerate and 

n 

ind/i|^ = ^indGi - indG. (3.23) 

i=l 

Proof. By Proposition 13.61 in the coordinates /ij = AAj G we have = xIq; 
where x is the following quadratic form on M"'^: 

n 

and the space @ C W^^ is defined by the condition 

i=i J 

Below we use the same notation x foi^ the corresponding bilinear form. Consider 
the fc-dimensional space 

E = {fi£W'' : in = GT^u, G M^ j = 1, . . . , n}. 
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Since G is nondegenerate, M"" = H © and moreover, xif^^O = for any € 
and £ E. Therefore, the spaces E and © are x-orthogonal and 

n 

ind Gi = indx = indxjc + ii^d x|e = ™d G + ind 

i=l 

□ 



Proposition 3.8. Y + Z = X andYnZ = T. 
This follows from 

Lemma 3.3. For any v € X there exists u = <l>v E y, unique modF, stic/i that 
u — V G Z. Explicitly, Ui = Vi + XjSiW^ , where the A^. satisfy (j3.16p and 

n 

= K^/^^ff^j^^/fl^^^.t^^+i), = = (3.24) 

1=1 

The map <1>: X — > y = y/F satisfies ^{Z) = and is the identity modF. 
Proposition 13.81 follows immediately. 

To prove Lemma 13.31 foi' given v € X we find u such that 

If{ui,Ui+i) = c" = const. 

n 

By Lemma [121 the X^i satisfy ()3.16p . If Xai is n-periodic in then AAqj = 0, 

i=i 

which gives ()3.24p . If (j3.24p holds, then equation (j3.16p determines n-periodic Xai 
modulo a constant independent of i. Thus, u is defined uniquely modulo F. 

Let X = X/T and \ei tty'- X ^ X be the corresponding canonical projection. 
Then there exists a linear map 11: X ^ X-^ such that the following diagram is 
commutative: 

X -X 

n 

X^ 

Corollary 3.2. The maps n|^ : Y — )• X-^ and *&|j^x • ~^ ^ ^^'^c mutually inverse 
isomorphisms. 

Proof. If IIv = for v G y, then v € Z and hence vGyn.Z^ = Fby Proposition l3.71 
The equation Y + Z = X implies that n|^ is surjective. We also have ^{Z) = 0. □ 
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The spaces Y,Z C X are orthogonal with respect to the bihnear form h on X. 
By Proposition 13.71 h is nondegenerate on Z and its index is given by (|3.23p . 

Prom the point of view of the Routh reduction it is natural to consider the space 

yO = {u G X : I,{ui, Ui+i) =0, i = l,...,n} CY. 

Define = (dai) G X* by 

dai = gaf5i-lBi-lWi_^ - ga/3iB*w^^^. (3.25) 

Note that {dai,wf) = 0. We have 

n 

{da,u) = ^gai3iIi{Ui,Ui+i), u G X. 

i=l 

Proposition 3.9. 

yo + Z = {v G X : (d«, v) =0, a = l,...,k}. (3.26) 
Proof. Let v G X. Then v G y° + Z, u = $(v) G Y^. By Lemma ESI 

n n 

= ^gal3iIi{Vi,Vi+l) = ^{dai,Vi) = (d„, v). 
i=l i=l 

□ 

Next we compute the restriction Let be the bihnear action form for 

the reduced system {E-^,A^), and let = ho where $ is the operator from 

Corollary [321 be the bilinear form on X"*- corresponding to h\y, see the commutative 
diagram 




Proposition 3.10. For any v G X , 

h^{^r,^r) = /i^(v, v) + 5„^c°c^ (3.27) 
where the coefficients c"(v) are defined by (j3.24p . 

This follows from a more general formula which we prove next. 
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Lemma 3.4. Let v £ X and u = ^*(v). Then 

n 

h{u,u) = /i(v,v) - go,pc°'c'^ + ^gapiIf{vi,Vi+i)I^{vi,Vi+i), (3.28) 

i=l 

where the coefficients c"(v) are defined by (j3.24p . 
Proof. By Proposition! 



h{u,u) = /i(v,v) + y^gf^AAgtAA^t- 

By (IHT^ and (ICTIl . 

n n 
^C/f'^AAaiAA/ji = ^5'a/3i(c" - If (t'i,t'i4.i))(c^ - /f (Uj, Uj+i)) 

i=l 

n 



2=1 i=l 

n 



□ 



Proo/ of PropositionlME We use (f3:28]l with v € X-^. Then 

If 77^+1 ) = -{BiVi,w°'_^-^), 

and so 

n n 

^9al3iIf'{Vi,Vi+i)I^{Vi,Vi+i) = ^gal3i{BiVi,Wij^i){BiVi,W^^-^). 
i=l 2=1 

Last we use (jS.lSp . □ 
3.7 The indices of h and h-^ 

In this subsection we discuss the relation between ind/i|y = ind/i|^ and ind/i^. 
Using the isomorphisms <I>: — )• Y and 11: y ^ X-*-, we compare instead the 
indices of and = ho ^ on X-^. As mentioned earher, h~^ and /i"*" coincide on 

We need some assumptions on the unit eigenspace of the linear Poincare map 
P: — 7- W. Suppose V = Ker(P— /) is isotropic. It is well known (see, for example, 
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|27| ) that the generahzed eigenspace = Ker(P — /)2m -g gyj^ipigctJc. Since V C N 
is isotropic, dimA^ > 2k. We consider the least degenerate case diniA^ = 2k. Then 

N = Ker(P - if = {v £ W : Pv - v £ V} 

is symplectic and V = Ker(P — /) is a Lagrangian subspace of N. Consider the 
bihnear form 

b{v, w) = uj {{P-I)v,w), v,w(£N. (3.29) 

A computation shows [27\ that b is symmetric. 

Let T, be the set of trajectories v of the variational system corresponding to the 
vectors v £ N. Then T C S. Trajectories in S \ T are not periodic. The projection 

n-^ = UT,C 

consists of periodic sequences. We have a natural map ^ = ^II: S — )• y. Set 
O = We will see that Jl-*" is orthogonal to = HY'^ with respect to /i"*", 

and O is orthogonal to Y with respect to h. Since and /i"*" coincide on X^, the 
difference of their indices is determined by their restrictions to the complement of 
Xq . Thus if O"*" © Xq = X-^ , the difference of the indices is determined by hJ | 
and /i^lj^x- However, in general X^ ^ ®Xq . To ensure this expansion we need 
an extra non-degeneracy condition O below. 

Take a basis ^ . . . m.V and the conjugate basis qi, . . . ,qk in a Lagrangian 
complement of V in A^. Then , . . . , qi, . . . , qk is a. symplectic basis in A^ and 

uj{w°',w'^) = 0, uj{qa,q/s) =0, uj{w" ,qfi) = 5'^, Pq^ = q^ + Sapw^ . (3.30) 

Combining (|3.29p and (j3.30p . we obtain 

Kla^qp) = Sap- 
ience the matrix S = (sa^) is symmetric. Define symmetric matrices A = (0^/3) 
and = (o^^): 

Oo;9 = Sa5H^'^Ssl3 — SajB, a^/S = Sap — Qap- (3.31) 

Below we need another non-degeneracy assumption. 
Condition C. The matrix A^ = [a^p) is nondegenerate. 

Definition 3.1. We say that a periodic trajectory is nondegenerate modV if the 
non- degeneracy conditions \^^[ and\C\ hold. 
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Theorem 3.2. Suppose that x is nondegenerate modV. Then 

ind/i - ind/ij^ - ind/i"^ = indA - indA^. (3.32) 
We prove Theorem 13.21 in ^ 13. 8[ 
Corollary 3.3. Suppose that x is nondegenerate mody. Then 

ind h = ind + ind + ind b + ind G (mod 2). 
Indeed, SA-^ = AG impHes 

sign det A = sign(det S det G det A^). 

Therefore, ind A — ind A-^ = ind b + ind G (mod 2). 
Proposition 13.71 implies 

n 

Corollary 3.4. ind /i = ind /i"*- + ind Gj + ind b (mod 2) . 

i=l 

Equation (j3.20p combined with Corollary 13.41 imply 
Corollary 3.5. cj(_i)md/i ^ ^±(_i)indh^+ind6_ 

3.8 The spaces Q and Q-^ 

Consider solutions 

w" = (wf), qa = (qai) 

of the variational system corresponding to the symplectic basis w", qa of the space 
A^. They form a basis in S. The solutions w"^ G F are n-periodic and satisfy (j3.12p . 
Equation (|3.30p implies that satisfy 

qa,i+n - Qai = Sa/BW^ . (3.33) 

Since the basis is symplectic, 

Ii'i<l/3i,qpi+i) = {BiW^, qpi+i) - {Biqpi, wf+i) = w(u'", g/3) = (5^, (3.34) 
{Bjqaj,q(3j+i) - {Bjqi3j,qaj+i) = uj{qa,q(3) = 0. (3.35) 

Let = Ilqa € il^ C X^. Then q^- = - X^pi'^i ^ where 

Xal3i+n ~ Xaf5i = ^o,l3, ^af^i = S-y/B i-l{Bi-lw]_^, qai) ■ (3.36) 
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For any a = 1, . . . ,k define Qq, G y by = ^'q^. Then 

where the coefficients Va/Si are chosen so that the Qa are periodic and 

IjimApj+i) = -{Bjqi3j,wf+i) + {Bjwf,q(3j+i) = (3.37) 
are independent of j. Then 

l^a0i+n — ^al3i = SafS, ^^afii = Soi-fH'^^ QSpi, = S'p — S ^.^K?" . (3.38) 

Since q^ is defined modulo F, we have q^ € 1" and q^ = ^q^ . 
We have the spaces 

= m: = span(q]^,... ,q^) C X^, = = = span(qi,. . . ,qfc) C Y. 

We define as the image of Y^ under the projection vrr : X ^ X = X/T. 

Theorem 3.3. Suppose that the non-degenaracy conditions [Aj [Bj and^hold. Then 

(a) Y = n@Y^ and X^ = ® X^; 

(b) the maps n|^ : — )• Q'^ and II\yo- Y^ Xq are linear isomorphisms; 

(c) h\yo = h^\^j_ ou\yo, /i|j.o = /i^l^x o n|yo; 

(d) for any u gY^ and a, f3 = 1, . . . , k, 

/i(qQ,,u) = 0, /i(q«,q/3) = ao/j; 

(e) for any E HY^ and a, /3 = 1, . . . ,k, 

/i^(q^,u^)=0, /i^(q^,q^)=a^^. 

The proof of Theorem 13.31 is contained in ^ lA.ll 

Now we prove Theorem l3.2i Recall that by Proposition [321 Y _L/i Z and Y ±^ Z. 
By Theorem 13.31 we have the /i-orthogonal expansion X = Z © e yo and the h^- 
orthogonal expansion X-^ = © X^ and /i|yo = ^^Ifiyo ° nj^o- Therefore, 

ind/i — ind/tj^ — ind/i^ = ind/i|^ — ind/i"'"|^x = ind(aa^) — ind(a;j^^). 
It remains to use the equations 

ind/i = ind/i, ind = ind ind = ind 
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3.9 Degeneracy for p-index form 

The connection between the indices for the original and the reduced systems is 
much simpler for the p-index form. We take complex p £ and perform the same 
computation for the corresponding Hermitian form hp on the complex space 

Xp = {u = {Uj)j^z ■■ Uj G Ej, Uj+n = PUj} 

of quasiperiodic sequences. We define Yp, Zp C Xp by the same formulae p.2ip : 

Yp = {ueXp: l\ui,U2) = ■■■= riun,Un+i)}, (3.39) 
Zp = {v £ Xp : Vi = XaiWf, Xai+n = pKj}- (3.40) 

The main difference is that for p ^ 1, w" ^ Xp. This implies, in particular, that 

Yp = Y^ = {u£ Xp-. Ij{uj,Uj+i) = 0, j = 1, . . . , n}. 

Proposition 3.11. The spaces Yp,Zp are h- orthogonal, that is, /i(u,v) = for all 
u £Yp and v G Zp. If /i(u, v) = for all v S Zp, then u € Yp. The restriction of 
hp to Zp is given by 

n n 

/i(v,v) = Y,9f^K^^hi = ^{GiAXi,A~Xi), 

1=1 i=l 
V = (AajW°), Xan+i = pXai- 



The proof is the same as for p = 1 (see Proposition 13. 6p . 

Suppose assumption \K\ holds Fl Then for p 7^ 1, 1"^ is the /i-orthogonal comple- 
ment to Zp and h\^ is nondegenerate. Indeed, 

/i|^(v,v) = (G,A,A), 

where Aj G C'^ and 

(GpA)j = Gj-iAAj-i — GjAAj, Aj+„. = pA?,. 

Lemma 3.5. 

n 

detGp = {-if p-^{p - if^WdeiG^. (3.41) 

i=l 



^We do not need assumption [Bl in this section. 
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This follows from Hill's formula (|2.18p applied to the DLS (F,Af) with the 
bilinear action form h\^. The corresponding Poincare map Pp has the matrix of the 
form 

7 

/ 



Pf 



Hence det(P — pi) = {p — 1)'^^. The operators Bi for the system {F,Ap) are equal 
to the Gi. 

Proposition 3.12. Let p e , p / 1. Then Xp = Yp® Zp and Ypr\Zp = {0}. 
Proof. Let u E Xp. If we want to find (A^i), Xai+n = P^ai-, such that 

(u, + A«iO G Yp, (3.42) 



then ()3.26p gives 



Thus 



1=1 



Aal = (/5-l) ^ y^ffaffi/f 
1=1 

Similarly we find Aa2 , • • • , Aq„ ■ □ 
Formula (|3.42p defines a projection ^p-. Xp Yp, ^pZp = 0, ^p\y = I. We 
have 

ind hp = ind hp\^ + ind /ip | y • 

The projection H : Xp — t- X^ gives an isomorphism of and the /3-index form /i^ 

for the reduced system (£'"'-, A"*-). Lemma 13.51 implies that the Hill p-determinants 
for the original and the reduced system are related by 

n 

det Hp = det det Gp = det H^{2 - p- p-^f JJ det Gi. 

i=l 

Hence 

n 

ind/ip|^ =^^indGj (mod 2). 

i=l 

Corollary 3.6. If p e , p^ 1, then 

n 

ind /ip = ind /ip + ind Gj (mod 2), null /ip = null /i^ . 
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4 Reversible version 
4.1 Reversible DLS 

Let S: M — > M be a smooth involution: S'^ = idjvf. We say that a DLS is S- 
reversible if S is a time reversing symmetry for L: for any x, y G M, 

L{Sx,Sy) = L{y,x). (4.1) 

Equivalently, the Lagrangian L is invariant under the involution S: M^, 
S{x,y) = {Sy,Sx). 

The simplest example is = id, that is, L(x, y) = L{y, x) (for example, a 
billiard system or a standard map). A non-trivial S appears in the system (j2.5p if 
the potential is even. Then S{x) = —x. An analogous possibility exists in billiards 
with some symmetry conditions. 

Proposition 4.1. Suppose that T: — ?• is generated by an S -reversible DLS 
and T{x,y) = (y, z). Then 

T{S{z),S{y)) = {S{y),S{x)). 
Hence T is conjugate to : 

T-^ oS = SoT. (4.2) 
The proof follows by differentiating the identity 

L{x,y) + L{y,z) = L{Sz,Sy) +L{Sy,Sx). 

If X = (xi) is a periodic orbit of a DLS, then x = (Sx-i) is also a periodic 
orbit. A periodic orbit is called reversible if x = x modulo translations. The group 
Z„ = Z/nZ acts on the set of n-periodic sequences {xi)i^z in M by translation 
(xj)jgz 1-^- {xi+k)i£Z, and we should identify periodic orbits obtained in such a way. 
Any n-periodic sequence is determined by (xi, . . . ,Xn) € and the translation 
group Zn acts on M" by cyclic permutations. Thus the set of periodic sequences is 
the quotient = M^/Z^. Define an involution R: — )• by R{x) = (Sx-i). 
Let ^+ = X € ^ : -Rx = x} be the set of fixed points of R. Thus x G if and 
only if S{xj-i) = Xi for some j € Z and all i (modn). 

Proposition 4.2. x € ^+ is a reversible periodic orbit if and only if x is a critical 
point of the functional = =2/ 1 . 

This is a well known property of functions invariant under an involution. Indeed, 
let X = Tx^. Then J = dR{x.) : X —t- X is an involution. Denote E± = € X : 
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= Then X = ® E_ and = T^^^. Since ^ is i?-invariant, we have 
d^(x) o J = (ij^(x). Thus, (ij/(x)^ = for all ^ G 

For any critical point x G let /i = d?£/ (x) be the Hessian, that is, the second 
differential of J^. This is a bilinear form on X. Then h\^^ = (i^,(S+(x) is the second 
differential of For ^ = ^+ + ^± G E±, we obtain 

Indeed, since /i is ./-invariant. 

If we represent /i by a linear operator H: X — >■ X*, then J*HJ = H, and so 
UE± = El. 

Let us introduce on M an S'-invariant Riemannian metric. It defines an R- 

invariant metric (•,•) on X. Then h{^,ri) = {HS^,r]), where H: X — t- X is a 
symmetric operator. The spaces E± are orthogonal with respect to the metric and 
HE± C E±. Denote H± = H\ej^. We obtain 

Proposition 4.3. H = H+ ® and det H = det i7+ det 

Reversible periodic trajectories x G ^+ are of 3 types r = 0, 1, 2 depending on 
the number of fixed points of S they contain. 

I^pe 0: n = 2k and x = (xi, . . . , Xj^, -Sxjt, . . . , Sxi). 

Type 1: n = 2A; — 1 and x = (xi, . . . , x^, Sx^, . . . , 8x2), where xi = Sxi. 

Type 2: n = 2k — 2 and x = (xi, . . . ,Xk,Sxk-i, ■ ■ ■ ,8x2), where Xi = Sxi and 
Xjfc = Sxf;. 

For all types x = XT-(y) G is determined by y = (xi, . . . ,Xk) G M^. Thus 
the action functional ^+ on ^+ gives a function on M*^: 

^^(y) = ^+(x^(y)) = ^+(x), r = 0, 1, 2. 

Denote 

fe-i 

i=l 

Lemma 4.1. x = y^riy) is a periodic orbit of type r = 0, 1, 2 if and only if y is a 

critical point of 

£/o{y) = 2^{y) + L{Sxi,xi) + L{xk, Sxk), y G M^ 

M(y) = 2^(y) + L(xfc,5xfc), yeNxM''-\ 

i^2(y) = 2^(y), y G iV X X iV, 

respectively. 
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The functional jz/ on ^ admits a similar representation. For example, consider 
the case of periodic orbits of type 0. A point x G M" can be written as 

x= {yi,...,yk,Szk,...,Szi). 

Then, since L{Szi^i, Szi) = L{zi,Zi^i), 

=(2/(x) = ^(y) + ^(z) + L{yk, Szk) + L{Szi, yi). 

For X a periodic orbit of type we have y = z = (xi, . . . , X}.). We write u G T^x^M'^ 
as 

u= JfcWfc,---,^itfi), v,w G TyM^, 

where Jj = dS{xi): T^^M — ^ TsxiM. Taking the second differential of we get 

h{u, u) = A;(v, v) + k{w, w) - {B^JkUUk, Vk) - {vi,BoJiWi), 

where 

k 

^(v, v) = ^{AiVi - Bi^iVi-i - B*Vi+i,Vi), vo = 0, v^+i = 0, 
1=1 

and Bk = B{xk,Sxk), Bo = B{Sxi,xi). Denote 

Ci = — -Bo Ji : El — >■ E*, Cfc = —B^Jk : E^ — >■ E^. 

Note that Ci = C* and = are symmetric. Thus 

h{u, u) = k{Y, v) + k{w, w) + {CkVk, Wk) + (Ci?^i, (4.3) 

Let us compute the corresponding bilinear forms h\^_^. For u G E± we have 
w = ibv, so u is determined by v: 

h±{u,u) = /i° (v,v) = 2A;(v, v) ± {CkVk,Vk) ± {Civi,vi). 

Similarly, for any r = 0, 1, 2 a vector u G E± is determined by v = (vi, . . . , vj^) G 

TyM. 

Lemma 4.2. For a reversible orbit of type r = 0, 1,2, h±{u,u) = /i^(v,v) /las t/ie 
form 

/t° (v,v) = 2/c(v,v) ± {CkVk,Vk) ± {Civi,vi), 
/i±(v, v) = 2A:(v, v) it {CkVk,Vk), 
4(v,v) = 2fe(v,v). 
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The domain of h'^ is V^, where 

Vl = {^r e TyM^ : JiVi = ±vi}, 

Vl = {^^e TyM^ : J^v^ = ±vx, JkVk = ±Vk}. 

Consider a periodic orbit of type 0. Then the domains of h\. and hP_ coincide 

and 

/i5i(v,v) - (v,v) = 2{CiVi,vi) + 2{CkVk,Vk). 

Corollary 4.1. Suppose that x G ^+ is a periodic orbit of type which is a 
nondegenerate local minimum point of s^^ . If the symmetric operators Ci and Ck 
are non-positive, then x. is a nondegenerate local minimum for . //x G is a 
nondegenerate local maximum point of and Ci, Ck are non-negative, then x is a 
nondegenerate local maximum for si . In both cases det H = det ff_ det > 0. 

For periodic orbits of type r = 1 or 2, the domains of h'!^ and /iL are different. 
When S = id, then C 1/^. 

Proposition 4.4. Suppose that S = idM, and letx be a reversible periodic trajectory 
of type T which is a nondegenerate minimum for ss/r- Then in each of the three cases 

(a) r = 2, 

(b) r = 1 and is non-positive, 

(c) r = and Ck, Ci are non-positive 

X is a nondegenerate minimum for s/. 

If y is a nondegenerate maximum for , then a similar statement holds provided 
that Ck and Ci arc non-negative rather than non-positive. 

Proposition 4.5. Let S = id. If si+ is a minimal periodic orbit of type 2, then hp 
is positive definite for all p £ S^ . Hence x is hyperbolic. 

Proof Take complex u G T^M"^^. Then 

hp{u,u) = A;(u+, u"^) + A;(u~,u~), 

where u+ = (ui, . . . , Uk) and u~ = {J2kU2k, ■ ■ ■ , Jk+iUk+i) are complex vectors from 
Since k is positive definite, hp is positive definite. □ 
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4.2 Some applications 

Corollary 4.2. Reversible geometric degeneracy det = of a reversible trajec- 
tory implies the dynamical degeneracy det(P — /) = 0. 

Next we give some statements on dynamical stability (in fact, instability) of 
reversible trajectories. 

Corollary 4.3. Let x be a reversible periodic trajectory such that y is a nonde- 
generate minimum point of Suppose also that cj(x)(— I)™' < 0. In the case 
type(x) = 1 we also need the condition that Ck is non-positive, and in the case 
type(x) = that Ci and are non-positive. Then x has a real multiplier > 1. In 
particular x is dynamically unstable. 

Proof. By Corollary 14.11 both h± are positive definite. Therefore det II± > 0. Now 
by (j2.8p and Proposition 14. 31 we have det(P — I) = —1. It remains to use Proposition 
[121 □ 

Consider, for example, the DLS generated by billiards in a domain in M™'^^ 
bounded by a hypersurface M. Suppose that the billiard hypersurface M is sym- 
metric with respect to a hyperplane in 1^™-+^, for definiteness passing through 0. 
Then the symmetry S = Se- W^'^^ — )■ M'"^-'^ is given by Se{x) = x — 2e{x, e), where 
e is the unit normal vector. 

Proposition 4.6. Let x, y S M , y = SeX, be a pair of symmetric points and let 
B = B{y,x). Then the operator C = —BSe- TyM — )■ T*M is symmetric and 
positive definite: {Cv,v) < for non-zero v G TyM. 



Indeed, by (I2i)]l . 

-{SeV,w) -\- {SeV,e){w,e) {v - e{v,e),w - e{w,e)) 

{Cv, w) = . . = . . , 

\x — y\ \x — y\ 

where we used SeV = v — 2e{v, e). 

Corollary 4.4. Let x be an Se-reversible periodic billiard trajectory of type t such 
that y is a nondegenerate minimum of the length functional . If m -\- t is odd, 
then X has a real multiplier greater than 1. 

Proof. Consider the case r = 0. By Proposition 14.61 the symmetric operators C\ 
and Ck are negative definite. Therefore by Proposition 14.41 x is a nondegenerate 
minimum of £^ . Note also that ct(x) = (—1)" = (—1)^ > 0. Now it remains to use 
Corollarv 12.61 The cases r = 1,2 are analogous. □ 

Any billiard is S'-reversible for S = id. Any reversible periodic trajectory x is of 
type 2. By Proposition EH a(x) = (-1)" > 0. 
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Corollary 4.5. Any id-reversible billiard trajectory which gives a nondegenerate 
minimum of the functional s^j^, is hyperbolic. 

This follows from Proposition I4.5i 



5 Hill's formula for a continuous Lagrangian system 

Consider a continuous Lagrangian system (M, ) with the configuration space M™ 
and smooth^ r-periodic Lagrangian ^(x, i, t) on TM x M. We assume that 5£ is 
strictly convex in velocity x € T^M. Then r-periodic trajectories are critical points 
of the action functional 




on the space of r-periodic W^'"^ curves 7: M — )• M. The goal of this section is to 
prove an analogue of Theorem 12.11 for continuous Lagrangian systems. 



5.1 Continuous Hill determinant 

The second variation of the functional at 7 € is a symmetric bilinear form 



h{^,r]) on the set X of r-periodic W^''^ vector fields £,{t) £ Et = T^u\M along 7. It 



is defined by 



do? 



a=0 

where 7^ : M ^ M is a smooth r-periodic variation of 7. Define a positive definite 
scalar product on Et by 

{v,w) = {B{t)v,w), B{t) = ^ii(7(t),7(t),t). 

Proposition 5.1. h can be uniquely represented in the form 

hiC,v) = f^[{Dm,Dri{t)) + {Uit)m,vit))) dt, (5.1) 

where U{t) = U*{t) is a symmetric linear operator and D is a covariant derivative, 
that is, a linear differential operator such that 

|(e,r/) = {Di,r1) + {i,Dl^), ^^{fC) = f^ + fD^ (5.2) 
for smooth vector fields ^{t),rj{t) £ Et and a scalar function f{t). 



^Actually, is enough. 
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Proof. Let V be any covariant derivative^ A standard computation shows that h 
can be written in the form 

hiC,v) = _^'(((Ve(t), V77(t)) + {Wit)at),Vv{t)) + {Vit)m,vit))) dt, 

where V{t), W{t): Et — )• Et are Hnear operators, and V{t) is symmetric with respect 
to the metric: V{t) = V*{t). 

By integration by parts h can be represented in the form (jS.ip . where 

Di = V^ + W -W\ U = V-W-W* = U*. 

Hence D is also a covariant derivative. □ 

Note that D and U are invariantly determined by h, that is, they are coordinate 
independent and do not change by a caUbration of the Lagrangian. 

Equations (|5.2|) imply that D is skew-symmetric relative to the scalar product 

iC,ri)2= r {m.ri{t))dt. 
Jo 

Therefore 

h{^, 7?) = {{-D^ + U)^, v)^ = (He, v)2, (5.3) 

where H = —D^ + 11 is the Hessian of £/ with respect to the L^-metric. 
The variational system of the periodic trajectory 7 has the form 

= U{t)m- (5.4) 

This is a linear Lagrangian system. We use the following definition. 

Definition 5.1. Let E = {Et}t<=R/TZ ^ smooth vector bundle. Suppose it is 
equipped with a metric ( • , • ) compatible with a covariant derivative D and a sym- 
metric linear operator U{t): Et — Et. Denote by {E,A) the linear Lagrangian 
system with the quadratic Lagrangian 

HtDO = \iDC,DO + ^m,0 (5.5) 



and Lagrange 's equations ([57 

* A covariant derivative is not uniquely defined: for an antisymmetric operator A{t), V + A{t) is 
also a covariant derivative. We use a covariant derivative because the derivative is undefined unless 
Et is ^-independent. 
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Trajectories ^{t), < t < t, of the system {E, A) are extremals of the quadratic 
action functional 

lh{^,0 = £ H^,DOdt (5.6) 

for variations with fixed ^(0),^(r). Thus h{$,,4>) = for any smooth (f){t) G Et such 
that <;^(0) = and 0(r) = 0. 

The system {E,A) is the linearization of (M, ^) at 7. In what follows we can 
forget about the non-linear Lagrangian system (M, =Sf) and work with the linear 
system {E,A). 

Let P : W ^ W he the linear Poincarc map of the trajectory 7. Since a solution 
^(t) of the variational system is uniquely determined by (^(0), Z)^(O)), W can be 
identified with Eq ® Eq. Then P is the monodromy operator of the variational 
system: 

P{m,Dm) = {^{r),DaT)). 
Define the VF^'^-scalar product on the Hilbert space X by 

((e, ^)) = {Di, Dr])2 + iC, v)2 = (BC, V)2, B = -L>2 + /. 

Then h{£,,fi) = {{H^,f])), where the self-adjoint operator H = B~^H is the Hessian 
of with respect to the M^^'^-scalar product. 

We have H = I+K, where K = {—D'^+I)~^{U —I) is compact, with eigenvalues 
Afe = 0{k~'^), A; = 1, 2, . . . , so that 

tr = ^ lAfel < 00. 
fe=i 

Thus, the Hill determinant 

00 

detH = l[{l + \k) 
k=l 

converges absolutely. 

Let Q: Eq ^ Eq he the monodromy operator of the equation of parallel trans- 
port: 

Q7?(0) = 7?(r), Dr]{t) = 0. 
Theorem 5.1. det(/ - P) = a(-l)"*/3deti?, where 

13 = e-"'^det2(e^/ - Q), a = det Q. (5.7) 
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Since Q is an orthogonal operator, /? > and a = ±1 depending on whether the 
bundle E is orientable, that is, if the trajectory 7 preserves or reverses orientation. 
If M is orientable, then a = 1 always. 

Theorem 15.11 follows from a more general result of the next subsection. 

Example 5.1. Suppose the Lagrangian system has one degree of freedom and the 
bundle E is trivial. Then Q = 1 and = — C + Since detP = 1, 

p-^ det{pl -P)= p + p-^ -2 + det(/ - P). 

If /9 is a multiplier. Theorem 15.11 gives 

detH= P + P~'-^ . (5.8) 

Set T = 27r and represent the operator H in the basis {e*"*}. If 

then 



Hence (ll.3p is the matrix of H, and (jS.Sh gives Hill's formula ()1.2p . 



5.2 Relation to Hill's formula for discrete Lagrangian systems 

A continuous Lagrangian system locally, near a periodic orbit 7, defines a discrete 
Lagrangian system. Take a partition = to < < ■ ■ ■ < = of [0, r] and let 
Xi = 7(ii). If the Ati = tj+i — ti are small enough, the points Xi, Xj+i are non- 
conjugate along 7. Then there is a neighbourhood Ui of in such that 
for each {x,y) € Ui there exists a unique trajectory Ux^y'- — > M close to 
7!, , and joining x and y. Define a discrete Lagrangian Lj on Ui by 

Li{x,y) = £/{ux,y). 
The discrete action functional A is 

n 

-4(y) = e f/j, yn+i = yi- 

i=l 
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Then x is a critical point of A, that is, the periodic orbit of the DLS corresponding 
to the periodic orbit 7 of the continuous Lagrangian system. It is easy to see that 

Bi = -di2Li{xi,Xi+i) = ^{B{ti)+0{Ati)). 

Thus Li is a discrete Lagrangian. 

The definition of the Hill determinant detH = det(B^^H) for a CLS is similar 
to (j2.8p for a DLS. Discretization of the operator H = + U corresponds to the 
operator H in (j2.8p . However, the operator B = —D^ + / does not correspond to 
the operator B in (j2.1ip . 

The choice of B is natural for DLS, but not so for CLS, where instead of / we 
could add almost anything. This is the reason for the strange coefficient /3 in (j5.7p . 
If we use an analogue of discrete B, then B^^H will be unbounded. 

5.3 Generalized Hill determinant 

For a given p G let Xp be the vector space of complex p-quasiperiodic locally 
W^''^ vector fields ^{t) G Et such that ^{t + r) = p6.{t)- Define a Hermitian p-index 
form i27j on Xp by ([5T]) : 

Ktf]) = £[{Dm,Dm) + {u(tm,m)) dt. (5.9) 

We also denote by X the complexification X = X'^, that is, the set of complex 
r-periodic T^^'^-vector fields along 7. For definiteness choose In p so that 

< Imlnp < 27r, 

and let fi = In p. Identifying X and Xp by the map 

X e^*C(t) G Xp, (5.10) 

we obtain a Hermitian form hp on X: 

hp{i,f)) = hie'^'^,^) = {{D + pI)^,{D + nI)rj)^ + {U^,fj)2 

= -{{D + ^iifi, n)2 + (t/e, n)2 = {{-{d + pif + u% r?)^ 

= (HpC,??)2, whereHp = -(Z) + ^/)2 + C/. 

(We used that p = —p and D is real and antisymmetric). Define the /o-Hessian 
operator Hp-. X ^ X hy hp{(,,f]) = {{Hp^,f])). Then 

Hp = B^Up = i-D^ + ly^ {-{D + pif + ^7). (5.11) 
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We henceforth assume that p 7^ may take any complex values. The operator 
Hp is self adjoint for \p\ = 1, but not in general. Although this is not a trace class 
operator: tr \Hp — I\ diverges for p 7^ 1, we can, following Poincare [2], define the 
generalized Hill determinant det Hp by means of the finite-dimensional approxima- 
tion 

det Hp = lim det Hf\ h\^^ = PNHpP^j : X^^) ^ x'^^\ (5.12) 

N-^OD 

where Pn is the orthogonal projection onto the finite-dimensional eigenspace X^^^ 
of the operator D corresponding to the eigenvalues in 

= {z. G A = a{D) : < N}. 

Theorem 5.2. The determinant ()5.12p converges and 

p-"" det{P - pi) = ai-l)"" 13 det Hp. (5.13) 

We present the proof of Theorem 15.21 in ^ 15.51 
For p = 1 we obtain Theorem 15.11 

5.4 Some applications 

Suppose that ^(x,i;) = (x,x)/2, where ( • , • ) is a Riemannian metric on M. The 
periodic orbit 7 is a closed geodesic. The quadratic Lagrangian of the variational 
system has the form (|5.5p . where D is the Levi-Civita covariant derivative along 7 
and U{t)^ = -R(C;7(0)C with R the curvature tensor. 

The variational system has a periodic solution ^{t) and a first integral {DS,, 7) = 

— (^,7) periodic in time. Hence P has two unit multipliers, and detH = det(I — 
at 

P) = 0. Let us present a reduced version of Hill's formula. More general results will 
be proved in the next subsection (see Corollarv lG.ip . 

Let = {u £ Et : {u,'y{t)) = 0}. If ^(t) G E^, then D^{t),U{t)^{t) G Et. 
Denote by if"*" : — ?> X-^ the restriction of H to the invariant subspace 

X^ = {^eX: m G E^}. 

Let P^ : W'^, = E^ x Eq, be the monodromy operator corresponding 

to solutions ^(t) G E^ of the variational system. Let : E^ Eq he the map of 
parallel transport along 7. Applying Theorem 15.11 to the linear Lagrangian system 
(i?-*-, A^), where A^ = A|^x, we obtain the following result [5]. 

Corollary 5.1. Hill's formula for the reduced system has the form 

det{P^ -I) = a{-ir-^P^detH^, = e-^^'-^^^det^Q^ - e^I). 
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Let us formulate another corollary to Theorem 15.21 For complex p let X-^ be 
the space of complex vector fields ^(t) G Et and let = ifp|^^. Then 

p-H-i) det(P^ - pi) = det . (5.14) 

Let us present a proof of ()5.14p from which will be generalized in ^ 16.11 Let 
Z = {ieX : i{t) = \{t)^{t)}. We write any ^ G X as i{t) = r]{t) + \{t)j{t), where 
7] £ X^. Then 

Kta = Kr],fj)+ r \X\^dt. (5.15) 
Jo 

Hence Hp = Hp | ^ © H-^ and det ffp = det Hp | ^ det i/^ , where det Hp | ^ is the Hill 

determinant for the system with the quadratic Lagrangian |Ap/2. The characteristic 
polynomial of the corresponding monodromy matrix is (p — 1)^. Thus by (j5.8p . 



e-(p-l) 
/o(e"-l) 



detF, 
But 

det(p/ -P) = {p- 1)2 det(p/ - P^), det(e^/ - Q) = (e^ - 1) det(e^/ - Q^), 

which implies ()5.14p . 

Next we discuss applications to stability of periodic trajectories, similar to the 
discrete case. For p € define the /?- index ([28], p7|) indp7 of a periodic trajectory 
7 as the index of the Hermitian form hp. Then ind7 = indi 7 is the Morse index of 
7. It equals the number of negative eigenvalues of the operator H. If p is not an 
eigenvalue of P, 

(_l)mdp7 = sign det i^p = a(-l)'" sign (p-™ det (/)/ - P)). 

The argument of the sign function is real for \p\ = 1 since the characteristic poly- 
nomial is reciprocal. 

The next result is proved in f5j. 

Corollary 5.2. Suppose the trajectory 7 is nondegenerate and a{—l)^~^^^'^'^ < 0. 
Then 7 has a real multiplier p > 1. 

Indeed, the characteristic polynomial F{p) = det {pi — P) satisfies F{1) < and 
F(+oo) = +00. Hence F has a real root p > 1. 

Corollary 15.21 is not true if 7 is degenerate. Suppose for example, that 7 is a 
closed geodesic o Then (j5.15p implies ind/ip = ind/ip , and so 

(_l)indp7 = sign det i/^ = (T(-l)'^-isign(/)i-'"det(p/-P^)). 



'The general degenerate case is discussed in §(5] 
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Corollary 5.3. Suppose the closed geodesic 7 is nondegenerate and 1)™^' > 
0. Then the characteristic polynomial F(p) = det{pl — P) has a real root p > 1. 
Therefore, 7 is exponentially unstable. 

This is proved in [5] and |29) using Hill's formula and also recently in [TT] using 
the theory of Maslov index [?]• In particular, nondegenerate closed geodesies of 
locally minimal length on an even-dimensional orientable manifold are exponentially 
unstable. Degenerate geodesies are linearly unstable, but in general instability will 
not be exponential and so has no relevance for applications to Lyapunov stability. 

Suppose m = 1 and let the 2r-periodic trajectory 7^ be 7 traversed twice. // 
7^ is nondegenerate (that is, ±1 are not multipliers), then 7 has hyperbolic (elliptic) 
type if and only if indj'^ is even (odd). 

Indeed, 7 and 7^ are simultaneously elliptic or hyperbolic. The multipliers of 
7^ are squares of the multipliers of 7. Hence, 7^ is hyperbolic if and only if its 
multipliers are real and positive, or, equivalently, signdet(/ — P^) = (—1)1+"^'^')' = 
— 1. Similarly, the ellipticity of j"^ is equivalent to (—1)1+™'^')' = 1. 

For the geodesic case, we obtain the following result of Poincare. Let 7 be a 
closed geodesic on a 2-dimensional Riemannian manifold. If 7^ is nondegenerate, 
then 7 has hyperbolic (elliptic) type if and only if ind7^ is even (odd). 

Suppose now p = —1. Then (—1)™^-!')' = fTsignF(— 1). Thus, if ind_i7 is 
odd, there exists a real multiplier p < —1. Note that the space X_i corresponds to 
antiperiodic variations such that ^(r) = — C(0). Since 2r-periodic vector fields are 
sums of r-periodic and r-antiperiodic, 

ind_i 7 = ind 7^ — ind 7. 

Thus, if ind 7 and ind 7^ are not even or odd simultaneously, then 7 is unstable. 

5.5 Proof of Theorem 15.21 

We follow [5J, see also [6j. The method goes back to Poincare's proof ([2], [3D]) of 
Hih's result p. _ 

The real skew-Hermitian operator D = D = —D* has compact resolvent {D + 
pl)~^. Its spectrum A = cr(L)) C iM coincides with the set of characteristic expo- 
nents of the equation Dr}{t) = of parallel transport. Thus 

K = {v: det(Q - e''^) = 0}. 

If € A, then —v and v -\- uj belong to A, where cj = ^Txijr. 

Let . . . , /3m be the roots of det(Q — pl^ = 0. Since \pj \ = 1, we may represent 
them as pj = e'^^'^ , where < IniVj < 27r/r. Then 

m 

A = [j{uj+ujZ). (5.16) 



50 



First suppose that fi ^ A. Then Hp = ST, where 

sifj.) = -i-D^ + iy\D + r(/i) = i-{D + iii)-^u. 

Since PnD = DPn, by ^J2\i we have 

det Hp = det S det T. 
The finite-dimensional approximation (j5.12p of the determinant 
detr(^) = /(//) = hm det(P^TP;^), 

N—i-oo 

converges absolutely for /u ^ A since 

tr |(Z) + ^/)-2[/| < oo. 

Hence / is a holomorphic function on C\ A having at points in A poles of multiplicity 
not greater than double the multiplicity of the corresponding points of the spectrum 
of D. 

The function / is periodic: f{i_i + uj) = /(/i). Indeed, if ^ G X'^, then e"^^^ G 

and 

{I-{D + iiI)-^U)e^'t = e"*(/ -{D + {ii + uj)I)-^U% 

so T(^) and T{fj, + uj) are similar. Thus f{n) = 4>{e^'^), where (j){p) is a meromorphic 
function having poles at the roots pi, . . . , pm oi det{pl — Q). The multiplicity of the 
pole is at most twice the multiplicity of the corresponding root. 

Hence there exists a polynomial g{p) of degree < 2m — 1 such that the functions 
(f){p) and g{p) det~^(p/ — Q) have the same principal parts of the Laurent expansion 
at each pole. Since (l){p) — )■ 1 as \p\ — >■ +oo, by Liouville's theorem, 

c^{p) = l + g{p)det-\pI-Q). (5.17) 

The determinant det S converges conditionally. By (j5.12p , 

dei{-{-D^ + I)-\D + plf)= hm TT ^^^±4^ 

!/GAjv 



/;/2 ,,2\2 ,,2 ,,2 

lin, (V)' TT (^^^) =(-l)Mim TT 



!/eAjv,«!/>o i/eA]v 
where A; is the multiplicity of zero in the spectrum of D. We have used that A = —A. 
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From (j5.16p it follows that the last product converges absolutely. Hence it is 
a holomorphic function of pi, ■ ■ ■ , Pm, P for pj ^ e^"^ and p 7^ 0. To compute the 
product, we will use Euler's formula (see, for example. 



n ' 



p \ cosh fiT — cosh UT 



{u + ojn)'^ J 1 — cosh i/r 



Suppose first that pj 7^ 1 and pi 7^ pj for i ^ j. Equivalently, Vj ^ and 
— Vj ^ ujIi for i ^ j. Then by (|5.16|) . 



1^2 



j=l neZ 



-r-J cosh pT — cosh l^jT -r-j- P + P ~ Pj ~ Pj 

^\ cosh T — cosh zy,r ■j--': + e""^ — p,- — 

^-rx ^''{p-P3? _ e"^-det^(p/-Q) 
p(e^ - Pjf p^ det^ (e^/ - Q) ' 

By continuity this holds for any pi, . . . , pm 7^ e^"^ and p 7^ 0. Hence 

p"^det2(e^/-Q) 

By and (ICTTI) . 

p'"detFp = (-l)V-ndet2(P^-Q)+5(p))- (5.18) 
Thus G'(p) = det Hp is a polynomial of degree 2m in p with leading coefficient 

We claim that the polynomials G{p) and F{p) = det(P — pi) have the same 
roots. It is sufficient to prove that if F{p) = 0, then G{p) = and the root has at 
least the same multiplicity. 

If F{p) = 0, there exists a non-zero r-periodic vector field ^ G X such that 
(—1)2 + C/)e^*^(t) = 0. Thus Hp^ = 0. Suppose first that p is not an eigenvalue of 
Q, that is, p ^ A. Then Hp = ST, where S is invertible and tr \T — I\ < 00. Hence 

= implies detT = (see, for example, [M]). Then det Hp = det 5 det T = 0, 
and so G(p) = 0. 
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If ^ € A, we can repeat the same argument replacing S and T, for instance, by 



5 = i-D^ + I)-^{-{D + filf + l), T = {-{D + filf + l) \-{D + filf + U). 

We have proved that 

G{p) = {-ifr'F{p). 

It remains to show that (—1)'^ = ct(— 1)™. Indeed, k is the dimension of the subspace 
on which the orthogonal operator Q is the identity, while a = (—1)", where n is the 
dimension of the subspace on which Q is a reflection. Since the dimension m — k — n 
of the complementary subspace is even, (j5.13p is proved. 

6 Degeneracy in Hill's formula 

In this section we consider the case when the periodic orbit 7 is degenerate, that 
is, the variational system has a non-zero r-periodic solution Equivalently, the 
linear Poincare map P has multiplier 1. Usually, this happens if the Lagrangian 
system has a time periodic first integral ^ which is nondegenerate on 7. Then, as 
proved by Poincare, the variational system has a non-zero periodic solution and a 
non-trivial linear time periodic first integral which is the linearization of ^ at 7. 
Here are two standard examples. 

1. Autonomous Lagrangian system. Then (M, ^) has the energy integral 

J^{x^x) = {p,x) — ^{x,x), p = ^±{x,x). 

The variational system of a periodic orbit 7 has a periodic solution = 7(t). A 
particular case is a closed geodesic in a Riemannian metric. 

2. A Lagrangian system with symmetry. Suppose the Lagrangian system 
(M, ^) admits a symmetry group ipg - M ^ M, s € M, preserving Let w(j;) = 

ipsix) be the corresponding symmetry field. Then 

s=0 

^{x,x,t) = {p,w{x)), p = J:fi{x,x,t), 



ds 



is the Noether first integral. The variational system of a periodic orbit 7 has a 
r-periodic solution ({t) = w(7(t)). Here is one concrete example. 
Planar 3-body problem. Here 

1 ^ mm 
i=i i^j ^^'1 
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Fix the centre of mass at the origin, so that 



M = Ix = (xi,X2,X3) G {R^f : ^m^Xj = 0, / xj L 



-1 

(Jxi, Jx2, Jxs), J = ( „ 1 . The Noether integral is the angular momentum 



Rotations of preserve and the corresponding symmetry field is w(x 

'0 -: 

1 



XyX^ — ^^{JXi^ miX{^ . 
i=l 

Since the system is autonomous, we have double degeneracy of any periodic orbit on 
which and ^ are independent: 1 is an eigenvalue of P with multiplicity at least 
4. Stationary periodic solutions (orbits of the symmetry group) have lower-order 
degeneracy. 

In several recent years many periodic solutions for the 3-body problem have been 
found by variational methods [13j. However, we see that the ordinary Hill formula 
is degenerate for them. In this section we put forward an approach to this problem. 

3. General Hamiltonian commutative symmetry. These examples are 
particular cases of Hamiltonian symmetries. Let us look at the system (M, ^) 
from the Hamiltonian point of view. Let be the Hamiltonian. Suppose that the 
system admits an algebra q of Hamiltonian symmetry fields v generated by integrals 
Let 7 be a r-periodic solution in the phase space and P the corresponding 
monodromy operator. For any v G ((t) = v(7(t)) is a periodic solution of the 
variational system. Therefore 

Pw = w, w = v(7(0)), V G 0. (6.1) 

If the system is autonomous, the Hamiltonian vector field of the system (M, ^) will 
be in g, and the corresponding eigenvector w is 7(0). 

In the present paper we consider only the case when the fc-dimensional algebra 
is commutative. Then the corresponding eigenspace V C Ker(P — /) is isotropic, 
and the multiplicity of eigenvalue 1 is at least 2k. Hamiltonian reduction makes it 
possible to remove this degeneracy, but then the reduced system loses the natural 
Lagrangian structure. 

The classical way to remove autonomous degeneracy is to pass from the Hamilton 
action functional to the Maupertuis action functional on the energy level [32j. The 
classical way to remove symmetry degeneracy in a Lagrangian system is the Routh 
method [32]. We briefly describe it here. 
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Suppose the system (M, ^) admits k commuting independent symmetry fields 
w^, . . . , w'^ on M: 

[w",w^]=0, a,P = l,...,k. 

The corresponding flows of symmetry tpf commute. Let G be the (local) commuta- 
tive group acting on M by X Vsi o • • • o i^sk (^)' ■s ^ Suppose that M = M/G 
is a smooth manifold and vr : M ^ M a smooth fibration. 

The Noether integrals = (p, w") give a vector integral ^ {x, x, t) G M'^. The 
Routh method reduces the Lagrangian system (M, ^) with fixed value ^ = c G M'^ 
of the Noether integral to a Lagrangian system (M, ^) on the reduced configuration 
space M. 

For c = the reduced Lagrangian is defined by0 

^ix, X, t) = min ^(x, x + So-w^'fx), t) , (6.2) 

provided the minimum exists, for example, ^ is superlinear in velocity. Since ^ 
depends only on x = Tr{x) and x = dTT{x)x, it can be regarded as a function on 
TM X M. 

For c ^ take closed G-invariant 1-forms vi, . . . on M such that fa(w^) = 
Then if we replace the Lagrangian by gauge-equivalent 

(x, X, t) = J£(x, X, t) — c"z^a(x), 

Lagrange's equations do not change, but the Noether integrals will be replaced by 
= — c". Hence the value c of the Noether integral is replaced by and so 
the Routh function can be defined by (|6.2|) . The following theorem folds (see 
[32]). 

Theorem 6.1 (Routh). Let x{t) be a trajectory of the system (M, ^) with ^ = c. 
Then x{t) = 7r(x(i)) is a trajectory of the system (M,^). Conversely, if x{t) is 
a trajectory of the system (AI,^), then there exists a trajectory x{t) of the system 
(M, ^) with ^ = c such that x(t) = 7r(x(t)). 

If 7 is a periodic orbit of the system (M, ^) and 7 the corresponding orbit of 
the system (M,^), then their variational systems are related by a linear version 
of Routh's method. In the next section we describe the Routh reduction for a 
linear Lagrangian system. It applies in a more general case, for example, when 
the Lagrangian system has non-Noether integrals. In particular, the linear Routh 
reduction includes the linearized Maupertuis reduction on an energy level. 

^"Recall that we assume summation in repeated Greek indices. 

^'^Such Ua exists globally if the fibration vr: M — >■ M is trivial. In general the first Chern class 
provides an obstruction. However Va always exists in a neighbourhood of a periodic orbit. 
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6.1 Routh reduction in a linear Lagrangian system 

If the linear Poincare map P : W W the periodic orbit 7 has eigenvalue 1, then 
to any eigenvector w = Pw there corresponds a non-zero r-periodic solution ((t) of 
the variational system {E,A). As proved by Poincare, the variational system has a 
linear r-periodic first integral 

Indeed, by ([O]) 

±ld^{t),Dm) = iC,D^() - i^,D\) = - (LUC) = 0. 

In fact, I({S,,D^) = u}{w,v) = Jw{v) is the value of the symplectic form on the 
vectors v,w £ W corresponding to ^, 

Suppose the Poincare map P has several eigenvectors corresponding to unit 
eigenvalue. Let V C Ker(P — I) be an isotropic subspace and let P C X be the 
corresponding vector space of periodic solutions of the variational system {E,A). 
Let w^, . . . be a basis in V and . . . G P the corresponding independent 
solutions. The variational system has first integrals 

/"(e, DC) = iC, DC) - iC, DC), a = 1, . . . , fc, 

in involution 

I''{C^,DC^) = iCDC^) - (C^, DC) = uj{w'', w^) = 0. (6.3) 

We write shortly / = {I^,. . . ,1^). 
Denote 

Ft = {at) : C e P} = 
To simplify the presentation we use the 

Condition A', dim Ft = k for all t. 

Equivalently, ^{t), . . . ,C{t) ^ are independent for all t. Thus the Gram 

^=(5"^), 9"^(t) = (r(t),C^(t)), (6.4) 

is nondegenerate for all t. 

In Appendix IA.3I we will show that this assumption is unnecessary. In fact, the 
set S = {t G M/rZ : dimF^ < k} is finite and the family {Ft)t^j^ can be extended to 
a smooth /c-dimensional vector bundle (-Ft)tgM/z- We will show that everything in 
this section works without the non-degeneracy assumption I A^[ 



--span{(:\t),...,C{t)}- 

following non-degeneracy assumption. 
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We describe Routh reduction for the linear system {E,A). The reduced config- 
uration spaces Ef = Et/Ft can be identified with 

E^ = {ueEt: {u, w) = for ah w G Ft} 

via the orthogonal projection H = Ht: Ft ^ E^. For a smooth field ^(t) € Ft 
denote D-^S,{t) = UtD{{t). Exphcitly, 

= e - Oapit C^)C, D^C = DC- QapiD^, C^)C, (6.5) 

where G^^ = {gaji) is the inverse of the Gram matrix G = (g"^). In Appendix lA. 31 
we show that 11 and are smooth also when the non-degeneracy assumption fails. 
Define the Routh Lagrangian A-L on E^ = {Ft) by 

A^(r/, D^r^) = ^ {D\ D^t]) + ^ {U^rj, rj), r/it), D^r]{t) G E^, 

where the symmetric operator U-^{t) : E^- E^- is given by 

{U^u,u) = {Uu,u) -3g^is{u,D^C){u,D^C^), u G E^ . 
Thus, = mj - 3C, where 

(C is independent on the choice of the basis). 

The bilinear action form of the system (£'-'-, A^) is 

]^h^{r,,'n)= K^{r,,D^rf)dt, v{t) e E^. (6.6) 

We have Routh's theorem for linear Lagrangian systems. 

Theorem 6.2. Let S,{t) £ Ef be a solution of the system (E, A) such that Z)^) = 
0. Then r]{t) = n^(t) G E^ is a solution of the system {F^,A-^). Conversely, if 
r]{t) is a solution of the system {E-^,K^), then there exists a solution ^{t) of the 
system (£', A), defined modF, such that I{^,D^) = and r](t) = n^(t). 

For the proof we need the following evident result. 

Lemma 6.1. Let 

ilit) - m G Fu m = V{t) + Xa{t)C{t). (6.7) 
Then ^ satisfies /f (C, F>^) = c" for all a = 1, . . . ,k if and only if 

K=gap{cP -L^{ri,Dr^)). (6.8) 
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Indeed, /"(e, D^) = I'^{r}, Dr]) + c/^'^A^^. 

Proof of Theorem \6.2l A vector field S,{t), < t < r, is a solution of {E,A) if and 
only if 

0) = r((^e, m + {ui, dt = 



for any smooth variation (j){t) G Et such that (/)(0) = (/>(r) = 0. 

Suppose I{^,DS,) = and let r] = 11^. We need to show that for every smooth 
variation (/}{t) G such that 0(0) = (/'(r) = we have 

h^irj, 0) = HiD^ D^cj)) + (;7^r?, 0)) dt = 0, 



where /i is the bilinear form (j6.6p corresponding to the Routh system. 
By Lemma Em C = ^ + -^aC"; where r/(t) € i?^ and 

Aa = -ga(sl^{r], Dt]) = 2g^^{r], DQ^). (6.9) 
Since (??, C") = (</>, C") = 0, I^M gives 

We obtain 

H^, ^) = rUDr] + XaC + KDC, Dct>) + {Urj, 0) + (A^f/^ , 0)) dt 
Jo 

= [{{Dr], D<P) + {Ut], <P) + (A„r + KDC, Dc/)) + {XaD^C, </>)) dt 
Jo 

= J^(^{D\D^<P) + {Ur],<p) + {g^^{r],DCK^,gse{<P,DC')C) 

+ ^((A«i?r, 'A) - 2(A,r, </>))) 

= r{{D^rj, D^(P) + {Ut], ct)) - 3g„;3(0, DC){rj, DQ^)) dt 
Jo 

Hence t? is a trajectory of {E^,A^). We skip the proof of the converse. □ 
Now we can write Hill's formula for the reduced linear Poincare map P: W ^ W. 

Let 

X^ = {7jeX: r?(t) G E^} 
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and let 

be the Hessian operator for the reduced system (£'-'-, A^). Let Q-^ : Eq — )■ Eq be the 
operator of parallel transport corresponding to D-^t] = and let o"-*- = det Q-^ = ±1. 
If assumption lAl holds, then a = because the bundle F is oriented. In general 
F can be non-oriented. 

Corollary 6.1. 

det(P-/) = a^(-l)™-*^^^ det H^, = e(™-*^)"det-2(e"/-Q^) > 0. (6.10) 

For the geodesic problem Q-^ = Q\j^±, = H\^^, P = P-^, a = a^, and we 
obtain Corollary 15. 11 

Note that in general ^ H\^^, except when D-^(^ = 0. The reason is that 
if 77 G X"*- is T-periodic, A in (j6.9p is not periodic in general, and so ^ ^ X. Hence 
the space X-^ of periodic rj{t) € E^- does not correspond to the space of periodic 
^(t) G Et such that /(^, D^) = 0. Thus is not the restriction of h to X^ as in 
the geodesic case. Hence we need to discuss the relation between and h. 

6.2 Elimination of degeneracy in the action functional 

As in ()3.2ip , define subspaces Y,ZcXhy 

Y = {^£X : /(e, DO = const}, (6.11) 
Z = {C£X : C{t) G Ft for aU t}. (6.12) 

If C^, . . . , G r are basis periodic solutions, then 

Z = |e(t) = Xa{t)C{t) : Xa{t + t) = A„(t), j\''^KXpdt < 00 

Lemma 6.2. For any ij G X there exists ^ = <I>(r/) G y, unique modF, such that 
— 7] € Z. Explicitly, r] = ^ — XaC"> where Xa satisfies ()6.8p with = c°'{r]) given 
by 

c^ = f,^P r g^slHv,Dv)dt, {k''^) = {g^p)-\ go.p= r Qapdt. (6.13) 
Jo Jo 

This follows from Lemma [6. II for periodic and r]. We have defined a projection 
$ : X — > y = Y/T which is identical on Y and <I> = on Z. We obtain 

Proposition 6.1. YnZ = r andY + Z = X. 
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Proposition 6.2. The spaces Y and Z are h-orthogonal. that is, ry) = for 

all G Z and rj . The restriction of h to Z has the form 



h{KCApC^)= g^^KXpdt. (6.14) 
Jo 

Proof. Take i^Z, i{t) = K{t)C{t). Then 

= HiKC, Drj) + {KDC. Dr,) + (A^Z^^C", r,)) dt 
Jo 

= (|(^, KDC) + A„r (r?, Drj)^ dt 

= {v,XaDC)\:+ r KI^{r^,Dr^)dt (6.15) 
Jo 

(we have used that Q satisfies the variational system). If r/ = A^C", then /"(r/, Drj) = 
g'^^Xp. Hence 

h{Xo.CApC^) = r g''f'K\pdt + {X^CApDC%. (6.16) 
Jo 

If T/ € y and the A^ are periodic, (j6.15p gives and (|6.16p gives (|6.14p . which proves 
Proposition I6.2i □ 

Let 

X = x/T, Y = y/r, z = z/T. 

Then Y ® Z = X . The bilinear form h defined a form h on X and Y Z, while 
is positive definite. 



Corollary 6.2. 



ind /i = ind , null /i = null + /c. 



Corollary 6.3. The projection 11: X ^ X^ defines an isomorphism 

n|p:y^x^, (n|^)-' = ^.|^^. (6.i7) 

Indeed, if = for i e Y, then ^ G Z and hence CGynZ = rby 
Proposition 16.11 Similarly, Y + Z = X implies that n(y) = X-^. 

Next we compute the restriction Let h^{r],r]) be the bilinear form for the 
reduced system {E-^ , ) . 
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Proposition 6.3. Let hJ = /i|^on ^ he the bilinear form on X-^ corresponding to 
h\y. Then for any rj € X-^ , 

h'^iij^T]) = h^{r],r])+g^pc''c^, (6.18) 
where the c°'{r]) = ^ = <5r/, are defined by (|6.13p . 

This follows from a more general formula. 

Lemma 6.3. Let ^(t) = r]{t) + Xa{t)C'^{t), < t < t, be as in Lemma [U7J[ Then 

Kt^) = H71,V) - r Qapil'^iv.Dn) - e){L^{r^,Dr^) - c") dt 
Jo 

+ {{XaC, XpDC^) + 2(A«r?, DC) + 2c" A„) j^. (6.19) 

Proof. We have 

h{i, = /i(r?, r?) + 2/i(r?, KC) + HKC, XpC^)- 
Now (f09D follows from (f6J5D . (fOBD . and ([63]). □ 

Proof of Proposition \6.3[ If r/ and A are periodic, then the boundary terms in (|6.19p 
vanish. By (j6.13p . 

9a/3 (c° - r (??, I?7?)) (c^ - I^iv, Drj)) dt 

= I gapL'^iv, Dr])I^{rj, Drj) dt - 2 I gapc'^I^iv, Dr,) dt + g^pc'^c^ 
Jo Jo 

ga^L'^iv, Dr,)I^{r], Dr,) dt - g^pc^c^ . 

Next we use ry = G X^ . Then I°'{r],Dr]) = -2{r],D^O, and so 

r 5a/3/"(??, Dr])L^{i^, Dr]) dt = 4 T g^p{r], D^C){v, D^C^) dt. 
Jo Jo 

Finally, 

h{v,v) = r {{D^r,,D^fi) + g^p{r],D^C){v,D^C'') + {Uri,fj)) dt 
Jo 

= h^{7], r/) + 4 r g^^irj, D^C){rj, D^C^) dt. (6.20) 
Jo 

It remains to substitute ([OO]) in ([6l9]) . □ 
Prom the point of view of Routh reduction it is natural to consider the space 

yO = G X : /(e, DO = 0} C F. (6.21) 
Indeed, by (^TM . h\yo = h^ o n|^o- 
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Proposition 6.4. 

yO + Z= |7?GX:^" 5a/3p^C^ r/) dt = o| . (6.22) 
Proof. We take tj ^ X. Then ^ = ^t] G provided that 

0=1 go,pI^{r],Dr,)dt = -2 [ g^p{D^C^ ,r,) dt. 
Jo Jo 

Here we have used that 

g^pl''{rj,Dr]) = ^^{g^^{rj,C'')) - 2g^p{ri,D^C^). (6.23) 

□ 

We see that + Z = X if and only if D-^C = 0. Equivalently, DZ C Z. Let 
X^ = = !^neX^: j'^ g^p{D^C^v)dt = o|. (6.24) 

Then X^ has codimension < A; in X-^. We have /i^ > and /i^ = on X^^. 
Since ind hJ = ind /i, 

ind < ind < ind h'^ + A;. 

Let = y/y'^. Since H: y — ?> X-^ is an isomorphism, dimfi < k. The integral 
/: y — > M*^ gives a map Q — ?> M'^. To compare the indices of and h^, in ^ 16.41 we 
construct a basis in ri, on which I is nondegenerate. 



6.3 Indices of h and h 

In this section we discuss the relation between ind/i = ind/i|y and ind/i^. Let 
P : W W he the Poincare map. As in ^ 13.71 we assume that 

N = Ker(P - I)"^ = {v e W : Pv - v e V} 

is symplectic and V = Ker(P — /) is a Lagrangian subspace in N. Let w^,. . . 
be a basis in V and qi, ■ ■ ■ ,qk a basis in a Lagrangian complement to V in A^. We 
define the matrix Sa/s by formula (|3.3Up and the matrices A = {aap) and A-^ = {aj^j^) 
by (|3.3ip , where x"-^ and g^i^ are the matrices in (|6.13|) . 

Definition 6.1. We say that 7 is nondegenerate modF if Assumption IA^ (p. [5^) 
/lo/ds and det 7^ 0. 
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Theorem 6.3. Suppose that 7 is nondegenerate modV . Then 

ind/i - ind/i-^ = indA - ind^^. (6.25) 

The formulation coincides with Theorem l3.2|, but the proof is different. We prove 
Theorem 16.31 in §[631 Since is positive definite, as in the proof of Coroharv 13.31 
we obtain 

Corollary 6.4. Suppose det A ^ 0. Then 
6.4 The spaces Q and 

Consider periodic solutions C^{t), ria{t) of the system {E,A) which correspond to 
t/;", Qa- Then the C"(*) periodic and satisfy (|6.3p . Equations (|3.30p imply 

r]a{t + T) -rja{t) = Sai3C^{t), {j^^, Dr]p) - {-qp, Dt]^) = Q (6.26) 

and {C,Dr)p)-{r^p,DC) = S1. 
For any a = 1, . . . , A; we put 

f]a = r]a- KjiQ^, 

where the coefficients Xap are chosen so that the r}^ are r-periodic and 

iCDf^p) - {f^p^DC) = 0% = const. 
Then the satisfy (|3.38p and 

Ao/3(t + 1') — Xapit) = Sap, = SaSK^'^Qep- 

We define t]^ = Hrja- Then 

Va = Va - Xaf^C^ , A^^(t + t) - A^^(t) = = {Va, C^)9Sp ■ 

Consider the spaces 

Q = span(r/i, ... ,fik) = ^il.^ CY, = span{r]i, . . . ,7]^) = UT. C X^. 

We also define Q and 1''' as the images of and Y'^ under the canonical projection 

Ur: X = x/r. 

Theorem 6.4. Suppose det^^ ^ 0. Then 

(a) Y = ne yo, x-^ = n^e hy^; 
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(b) the maps n|^ : J7 ^ $7 and II\yo' Xq are linear isomorphisms; 

(c) hlyo = h^\x^ °n|yo and h\yo = h-^\x^j- °'^\yo' 

(d) for any ^ gY^ and a, (3 = 1, . . . , k, 

(e) for any S HY^ and a, (3 = 1, . . . ,k 

h^{v^,e) = 0, h^{rj^,rjj) = a^f,. 

The proof of Theorem 16.41 is contained in ^ IA.2i 
Corollary 6.5. In the basis fji, . . . ,fiis 

h\^-h^oU\^ = {sa,sK^'Se0-2Sal3+9ap)= SKS -2S + G. 

This quadratic form is positive definite. 

Indeed, let Q be the square root of K, that is, the positive definite symmetric 
matrix such that = K. Then Q^^ = G and 

SKS -2S + G = RR*, R={S-G)Q. 

This matrix is positive definite because R is nondegenerate. 

Now we prove Theorem 16.31 Recall that by Proposition 16.21 and Theorem 16 . 41 we 
have the /i-orthogonal expansion X = Z (B^(B Y^ and the /i-'--orthogonal expansion 
= Q-^ © HY^. The form h\z is positive definite and /i|yo = ^"""Ix-l ° nj^o- 
Therefore, 

ind/i — ind/i"*" = ind/i|j^ — ind/i"'"|^x = ind A — ind^"*". 
6.5 Example: autonomous systems 

Suppose the Lagrangian system is autonomous, so the variational system of a peri- 
odic trajectory 7 has a periodic solution = 7(t). If 7 is nondegenerate in the 
autonomous sense (only two unit multipliers) then, as proved by Poincare, there 
exists a famil3@ 7q, of r(Q)-periodic orbits such that 70 = 7 and t(0) = r (j33|. 

[53]). Let E{a) = Jif\ and A{a) = / {p,dx) be the energy and Maupertuis 

action of 7^. 

^^It seems more natural to parametrize the family by the period r. However, this is not always 
possible because it may happen that r'(a) — 0. 
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Lemma 6.4. Suppose that dE/dr ^ 0. Then 

dV 

(_l)mdfe ^ _ sign(T'{a)E'{a)) = - sign ^ . 

' dr 

Proof. The union of trajectories of 7a in the phase space TM = T*M is a symplec- 
tic cyhnder S. Restricting the Hamiltonian system to S we obtain an integrable 
Hamiltonian system with one degree of freedom and Hamiltonian H{'d,I) = E{I), 
where ■(? € M/Z, / € M. Then we can assume a = I, "fait) = {E'{a)t,a). Then 

r(a) = l/u{a), where the frequency z/(a) is E'{a). We have w = and v = . 
The monodromy matrix of 7^ is P = Pq, = ( ^ ^ (^] . Thus Pv = ('^ {'^] , and so 



^0 1 y \ 1 ^ 

s = uj{Pv,v) = u' /v = -T-'^ dr/dE. □ 

As usual, we denote ind7 = ind/i(7). 

Proposition 6.5. Let a periodic trajectory 7 have exactly 2 unit multipliers. Sup- 
pose that a{—l)^~^^^^"'dE/dT < 0. Then 7 has a real multiplier p > 1. 

Proof. Since a = o""*", by Corollary 16.41 and Lemma [67 



(_^yndh^ = -sign(^^^(-ir'lT. 
The dimension of the reduced system is m-*- = m — 1. Hence 

and by Corollary 15.21 applied to the reduced Hill formula (I6.10p . there exists a 
multiplier p > 1. □ 

Example 6.1. Suppose a particle in moves under the potential field with ho- 
mogeneous potential energy 

V{Xx) = X''V{x), A>0, k{k- 2)^0. 

Suppose 7 is a r-periodic solution with energy E. Then jx{t) = \-i{\^'^-H) IS a 
periodic solution with period t(A) = A^^^/^r and energy E{X) = X^E. Hence 

dE{X) 2k /^(A)y^'+2)/{^-2) 

-M>^ = k^2[-^) ('-'^^ 

(-ir^'' = sign^. 



Thus by Lemma 167 
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Consider the problem of the motion of a particle in M"^ in the force field generated 
by a homogeneous potential of degree k, where k{k—2) ^ 0. Equations cr = 1, ()6.27p . 
and Proposition 16.51 immediately imply 

Proposition 6.6. Let a periodic trajectory 7 have exactly 2 unit multipliers. Sup- 
pose that ( — — 2)/k < 0. Then 7 has a real multiplier p > 1. 

6.6 Degeneracy in the p-index form 

We have seen that the relation between ind/i and ind/i"*" is not evident. This sim- 
plifies drastically for p / 1. Let Xp, p ^ S^, he the set of complex p-quasiperiodic 
vector fields. Similarly to (|6.1ip . (|6.12p define subspaces Yp,Zp C Xp-. 

Yp = {^GXp: /(C, DO = const}, Zp = e Xp : ^t) G Ft}. 

It is easy to see that for p 7^ 1 and ^ G Xp, /(^, DS,) = c implies c = 0. Thus, 

Yp = Y^ = {CeXp:I{^,DO^0}. 

Proposition 6.7. For p ^ 1 we have YpH Zp = {0} and Xp = Yp® Zp. 

Proof. We will define a projection $p: Xp — > Yp along Zp. Take rj G Xp and look for 
Xait) such that ^ = ry+A^C" G ^p- Then by ([62]), A„ = fa, where fait+r) = pfa(t). 
Hence fa{t) = e'^^ba{t), where p = Inp and ba{t) is a r-periodic function: 

i — ^ 7- 

We obtain a unique solution \a(t) such that Aa(t + r) = pXa{t): 

The denominator is non-zero if p 7^ 1. □ 

Proposition 6.8. The spaces Yp, Zp are h-orthogonal: 

HC, V) = for all Zp, r]£ Yp. 
The restriction of h to Zp is positive definite for p ^ 1: 

h{X^C,\pC^)= r g'^^XoXpdt. (6.29) 
Jo 
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Proof. Take ^{t) = A„(t)C°(t) G Zp and r] € Yp. Then by (l6T5]l . 

/i(A„r,r?) = {X^DC,f])\l = {\p\^ - l){X^{0)DC{0),Dr]{0)) = 0. (6.30) 

The proof of ()6.29p is similar. □ 
Next we compute the restriction of h to Yp. Let 

= {r]^Xp: G E^}, 

and let 11: Xp — )• be the projection (j6.5p . Since Xp = l^p = © Zp, 
n|y : Yp ^ X^ is an isomorphism and its inverse is $p|^x. 

Proposition 6.9. For p ^ 1 the bilinear form hj = h o (nj^ ) ^ on X^ is equal 
to the Routh form hj; . 

This follows from Lemma |6. II (for complex vector fields) since c = and A(t) = 
pA(0), 7?(r) =p7?(0), \p\ = l. 

Corollary 6.6. For p ^ 1 the p-index of the system (E, A) equals the p-index of the 
Routh system {E-^,A^). 

Proposition 16.91 is not true for p = 1. Then the relation between indices is more 
complicated, as we saw before. 

7 Reversible case 

Suppose the Lagrangian system (M, ^) is reversible: there is an involution : M — ?• 
M which is a time reversing symmetry for ^: 

^{S{x), dS{x)±, t) = if (x, -X, -t). 

Let T = 2T. Then for any r-periodic curve 7 G ^2, 

^ij) = l'^^^ij{t),^it),t)dt = ^ij), 

where ^{t) = Sj{—t). Thus, the involution R: $7^0, 

R{im = j{t) = s{7(.-t)), 

preserves A r-periodic orbit 7 is called reversible if -R(7) = 7. Then 
S7(-t)=7(t), S^{T-t) = ^{T + t). 
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Hence 7(0) and j{T) belong to the set N of fixed points of S. It is easy to see that 
7 is a reversible periodic orbit if and only if 7+ = 7|p^j is a critical point of the 
action functional ^ 

^+(i/)= / ^{v{t),v{t),t)dt 
Jo 

on the set f2+ of curves v. [0,T] ^ M with end-points in N. 

Let X = TjQ, be the set of vector fields along 7 and J = dR{^) : X ^ X. Then 

{Jm) = J-t^i-t), Jt = dS{j{t)) : Et ^ E_t. 

Since R preserves s/, the involution J preserves the Hessian bilinear form: 

h{J^,Jri) = h{^,v). 

Since the operators D and U are intrinsically associated with h, 

J* = J, DJ = -JD, UJ = JU. 

Let X± = e X : = ±^}. Then X = X+®X^ and any ^ G X is represented 
as ^ = ^+ + ^_ , where ^± G X± . Then 

Since DJ = — JD, we have D^: X± X±. Hence the Hessian operator H = 
(-1)2 + J)-i(-£)2 _^ jj-j commutes with J, and so H : X± ^ X±. Denote H± = 

Proposition 7.1. H = H+® and det H = det H+ del iJ_ . 

Next we give more explicit formulae for h±. Any ^ G X± is determined by the 
restriction 

C|[o,T] ^Y± = {r^eY: r?(0) G 7?(r) G F = W'^^i[0,T], E), 

where 

£J± = {t, G ^0 : Jot, = ±v}, E^ = {v€Et: Jtv = ±v}. 
Thus we have the orthogonal decompositions 

Eo = E+® E^, Et = E+® E-. 

For r] eY± the corresponding ^ G X± is given by 

^l[o,r] ^ ^l[-T,o] ^ 
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Thus, 

Ki, = KV: r]) + h{j7], Jr,) = 2K{7], 77), 

where 

K{V,V)= {{Dr],Dr]) + {Ur],r]))dt 
Jo 

is the same form h, but considered on Y. Let K± = K\yj_- Then = d^£^(7-|_) 
is the second variation of the functional . 

Let us consider the case 5 = id. Then Eq = Eq, E^ = Et, Eq = 0, Ei^ = 0. 
Thus, Y+ = Y and 

Y^=Yo = {rjeY : r/(0) = 0, r/(T) = 0}. 

Corollary 7.1. Let m be odd and S = id. // 7+ is a nondegenerate minimum of 
A+, then the corresponding reversible periodic orbit 7 has a real multiplier > 1. 

Indeed, 7 preserves orientation, so o" > 0. The Hessian is positive definite, 
and hence the same is true for K- = 

An analogue of Corollary [7T] is true also for m even. 

Proposition 7.2. If S = id and 7+ is a nondegenerate minimum of then hp 
is positive definite for \p\ = 1. Hence there are no multipliers on the unit circle. 

Proof. For complex ^ G X, set 

7?(t) = e'^^^it) = u{t) + iv{t), u{t),v{t) G Et. 

Then 

hp{i,i) = h{r], ff) = h{u, u) + h{v, v). 

Let us show that h{u,u) is positive definite on W'^''^{[—T,T],E). Indeed, u+ = 
G Y and n_ = J(n|j G Y. Thus 

h{u, u) = K{u^,u^) + K{u^,u ) > 0, li 7^ 0. 

By Hill's formula det(P — pi) ^ for p G S^, and so P has no multipliers on S^. 
For another proof see [35] . □ 

Consider again the case of a general involution S. Then n y_ =Yq. Let Y^ 
be the i^T-t-complement of Yq in Y±, that is, the set of G 1± such that K±{ri, (") = 
for all ( gYq. By integration by parts, 

y± = y± n y^, = {?? g y : (-d^ + u)r] = o}. 
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The restriction of K± to equals 

Ki{r,,r^) = {Dr,,r,)\l, r/ G Y^. (7.1) 
Let Kq = K\y^ and = K\y^. Then we have 

It foUows that if deti^o / 0, then signdet/7 = (_i)md/r|+indi^f ^ (_i)mdir^_ 
Proposition 7.3. If the time moments and T are non-conjugate, then 

(_l)ind7 ^ signdeti? = {-if""^^^. 

If the time moments and T are non-conjugate, dimY^ = 2n and dimyj- = 
2(m — n), where n = dim A^. 

The quadratic form has a simple meaning, the Hessian of the discrete La- 
grangian (Hamilton action function) defined locally as 

L{x,y) = 

where : [0, T] — )• M is a trajectory joining x and y. 

A Appendix 

A.l Proof of Theorem D 

(a) By (|3.37p . for any constant vector A'^ we have: 

/"(A%„A''<7^,+i) = c^A'^. 

By condition [Cl the matrix Ca is nondegenerate. Therefore the equation 

/°(A%i,A%,+i) = r" 
with respect to A'^ is solvable for any constant vector r". This implies the first 



statement in (a) 



To prove the second statement in (a) , we show that dim = k and 







0. In view of equation (j3.26p . the last two conditions are equivalent to the non- 
degeneracy of the matrix 
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By (|3.36p we have 

i 

= ^ 9a5i {Sp + ^peidf - '^jiei+ldi^) 

i 

= Qap - ^(A^£j_,.i - A^^j) = gafi - Sap = -a^p- 

i 

|(b)| Since the map n|^: 1" — )• X-^ is an isomorphism, n|^: $7 Q-^ and 
n|^>o : 1"*^ — > Xq are isomorphisms. 



(c) By Proposition 16.31 for any v € X 

/i^(v,v) =/i|^on-i(v,v)+5„^c"c^. 

If V € Xq-, we have v = IIu = (uj — Xaiwf), u E y'^. Then by Lemma [67 
c" = If (uj, Uj+i) = 0. This imphes the first equation in (c) To prove the second it 
is sufficient to note that h = ho tty- and 11 = vr o Ilr . 



(d) By (mni), 



i 
i 
i 

= Sa-yK'^^{B*gspiw'^^^ - Bi^igsp i-lW^_i, Ui) = -Sa^K'^^{dsi,Ui), 
i i 

where the dsi are defined by (j3.25p . Proposition 13.91 imphes the first assertion in 



(d) Analogously, 

^(qa,q/3) = SayK^^ "^{B* gseiwl^i - Bi_igsei-iwl_i,qj3i) 

i 

= Sa-yK'^^ ^ g5si{{Biqpi,wl^i) - {BiW^ , qp i+l)) 



This imphes the second assertion in|(d) 
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(e)]We get 



, w. 



i+l/ 



i 

i 

i 

= ^ga5i{Biui,wfj^^). 

i 

Since uj- = Ui - gefSi-.i{Bi^iwl_-^,Ui)v^ , we obtain 



This implies the first assertion in (e) Analogously, 



^^(q^.q^) = '^{9a5i{Biqi3i,wf^i) - ga5i-l{Bi-lwf_,i,qi3i)) 
i 

= '^{-gai3i + ga5i{Biwf,q/3i+i) - gaS i-l{Bi-lwf_i, qi^i)) 



-9aB + SaB- 



This implies the second assertion in (ej 



A. 2 Proof of Theorem 16.41 

(a) For any constant vector we have I'^{X'^fji3,DXf^rjfs) = c^X^. therefore for any 
constant vector the coefficients A'^ can be chosen so that I"'{X^fip, DX^fii3) = r". 



This implies the first equation in (a) 



To prove the second equation in (a) , we show that dim = k and Q^CiXq- = 0. 
In view of equation (|6.24p , it is sufficient to check that the matrix 

ga5{D^C\4)dt 
is nondegenerate. We have: 

Jo Jo 

{gas{DC',Vp) - 9as{Vp,C)9ue{DC',C)) dt. 





72 



Using the equation 



= I {geug"') = {geuC + 2g,,DC, C') , (A.i) 



we continue: 

I rfd 



/o [di^'^"^^'"^'^'^^^ - dap] dt= - {Sal3-9a(3) = -^CL^fj- 



|(b)| The maps n|^: — > and n|^o: y*^ ^ Xq are isomorphisms because 
n : y — 7- X-^ is an isomorphism. 



(c) By Proposition 16.31 for any rj G X 

h^{v,v) = h\y on'\r,,ri) + g^^c'^c^ . 

If 7] € Xq, we have ry = 11^ = ^ — A^C", ^ € Y^. Then by Lemma 16.11 c" = 
= 0. This imphes the first equation in (c) To prove the second it is 
sufficient to note that h = h o t^y and 11 = 11 o vrp. 
1(d) I Integrating by parts we get 

KfjcO = [ i-D^a + Ufj^,0 dt = s^s^^' [ {geuC + 2ge.DC, dt 
Jo Jo 

= sas^'' j^{gsu{C,O)dt = 0. 

Analogously, 

HVaiVp) = SaSl^^^ / {geyC + '^QeuDC : rip) dt. 

Jo 

We define fia/s so that fia/s = c^gsfs- Then 

(r , D{fip - fipsC')) - - f^psC', DC) = c1- fipsg'" = 0. 

Therefore, 

£(9suC + 2g,,DC,fip)dt = ^^{g,,(C,fip - ^psC')) dt 

Hepdt = -gescp, 

where we have used (|A.ip . Finally, h{fia,fip) = —SaSK^'^geuC'^ = —Sasc^ = (lap- 



(e) Integrating by parts we get 



iVa ' ) = {{{-DUD + U)vLC^)- Sgse {Va , DCs){^^ ,DC,))dt. 
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Direct, but lengthy computation gives 

Here we have used (|A.ip . Analogously, 
Jo Jo 

= j {^QaP + ^^{ga^iC^ dt = -gap + Sa/3- 

A. 3 Degenerate case 

In this subsection we consider the case when the nondegeneracy assumption IA1 on 
p. [56] fails, that is, rank(Cni), • • • , C''(*)) drops on S C M/rZ. We will see that the 
Routh reduction of the system {E,A) to (£'-'-, A^) and other results on elimination 
of degeneracy hold with minor modifications of the proofs. Note that for DLS 
condition [Aj which is similar to IA^ is probably necessary. 

Lemma A.l. The family {Ft)t^s can be extended to a smooth vector bundle F = 
{Ft)t(z^/TZ- Thus the orthogonal complement E-^ = (E^-) is a smooth vector bundle. 
The operator D-^ on E-^ and the reduced Lagrangian A^ defined for t ^ E can be 
smoothly extended to t ^T,. 

Proof. Suppose that E S, and let 

rank(CnO),...,C'(0)) =A:-/. 

Without loss of generality we may assume that 

C^(0) = • • • = C\0) = 0, rank (C'+^O), • • • , C'(0)) = k - I. 

Then rank (Z)(^^(0), . . . , DC'^(O)) = or else the solutions C . . . of the variational 
system are dependent. Since 

r (C^ DC^) = iCDC^) - ic", DC) = 0, 

we have 

{DC{0)X^{0)) =0, a = l,...,l, P = l + l,...,k. 

Thus 

rank(DCHO), • • • , I^C'W, C'+^(0), . . . , C'(0)) = k. 
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Since 

D{C - tDC) = -tD'^C = -tUC = 0{t^), 

we have 

Q''{t) = tDC{t) + 0{t^), a = l,...,L 

Thus, the space 

Ft = spanQ DC\t), . . . , J DC\t), C'+\t), C\t)j , t ^ 0, 
has a hmit 

Fo = span(Z)C'(0), . . . , DC'{0), C'+\0), (^=(0)) 

as t ^ 0, and Ft is smooth at t = 0. The first statement is proved. 

Since E^- is smooth at t = 0, also Ht: Ft — >■ E^- is smooth, and hence the 
operator D-^ is smooth. Finally we need to check that the term 

(Cr/,r/)=5a/3(u,DV)(«,^^C'') 

in A"*- is smooth at t = 0. 

Denote j"'^ = (L>C"(0), £'C'^(0)) , a, /3 = 1, . . . ,1. Then (j"/^) is a nondegenerate 
matrix and 

{Dc,Dc'') = r(' + o{t'). 

Thus 

gaf^ = (C", C'') = t'{DC, DCf") + 0{t^) = t'r^ + 0{t^), a, /3 = 1, . . . , L 

The matrix {g°'^{0)), q, /3 = / + 1, . . . , A;, is nondegenerate, while one can show that 

g^P = 0{t^), a = l,...,l, (3 = l + l,...,k. 

Thus for the inverse matrix (^o^) we obtain 

Gap = t~^{jai3 + 0{t^)), a,P = l,...,l. 

The block (ga^), Oi,/3 = I + I, ■ ■ ■ ,k, is smooth and nondegenerate and the block 
{gap)i a = 1, • • • , /3 = / + 1, . . . , A:, is smooth. 

Since D^C,°^{t) = 0{t^) for a = 1, we obtain that C is smooth at t = 0. 
Thus the reduced Lagrangian A"*- is smooth on E-^. □ 

In fact, everything we have done in §[6] holds in the singular case. For example, 
let us check that the projection X along Z is well defined and smooth. As 
in the nondegenerate case, we have 

Z = {ieX: m e Ft} = = Xa{t)C{t) : j^'^^KXpdt < ooj. 
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but now Aq, may be singular for t € S. 

Take r] £ X and look for Xa{t) such that ^ = rj + AqC" € Y. Then by (j6.8p and 

(IS231), 

The last term is smooth for t € S. Suppose again that € S. Then we obtain 
1 ' 1 ' 

K{t) = --^ia/jc'^ + -X]-?c./3(^,C^) + smooth terms, a = l,...,l. 

P=l 13=1 

It follows that \a{t)C°'{t) is smooth, so ^{t) is smooth. 
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